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SYMBOLS 


A, A 
B 
b 
E 


Ec 



k 

L 

m 


scale constants appearing in the interaction potential 
rotational energy constant 

impact parameter, or miss distance at infinite separation 
energy 

. mu)^L^ 

characteristic energy, — ^ — 

0.6581x10"^^ (eV-sec) Planck’s constant 
Boltzmann constant, ^l 6^05 4 

characteristic range of the interaction potential 
reduced mass of the collision 


P,P(Av,AZ) probability that a single collision with energy E and miss 

distance b will cause a transition {y ,1) ^ (v+Av,Z+AZ) 

r(t) distance between the oncoming particle and the diatomic 

molecule 


s 

S,S(Av, Al) 

T 

t 

U 

u 



symmetry factor 

cross section for the transition whose probability is 
P(Av,Al) (see eq. (4)) 

temperature 

time 

interaction or perturbation potential 

relative velocity of collision partners at infinite 
separation 

average relative velocity at temperature T 


X = 


mu^ 

2Ff 


dimensionless collision energy 


total rate coefficient for vibrational transition v^v+Av 


iii 


a(Av) 



a(Av,Al) rate coefficient for the transition (v,Z) ^ (v+Av,Z-+Al) 

3(Av,AZ) perturbation matrix overlap integral (see eq. (42)) 



y 


reduced mass of the molecule 


P 

Pe 

a 


O3,C0 (Av, AZ) 


separation distance of molecular nuclei 
equilibrium separation distance 

distance of closest approach, or turning point (see fig. 1) 

vibrational relaxation time (see eq. (6)) 

wave function of the rotating vibrating diatomic molecule 

time independent eigenfunction 

circular frequency 
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THREE-DIMENSIONAL COLLISION INDUCED VIBRATIONAL 


TRANSITIONS IN HOMOGENEOUS DIATOMIC 
MOLECULES 

C. Frederick Hansen and Walter E. Pearson 
Ames Research Center 


SUMMARY 


The one-dimensional, semiclassical theory of vibrational transitions in 
diatomic molecules is extended to three dimensions. Simple exponential inter- 
action potentials are assumed and are spherically averaged to determine the 
collision trajectory that defines the perturbation. For use over the range of 
temperatures where small perturbation theory applies, fully analytic approxi- 
mations are derived for the cross sections, the rate coefficients, and the 
relaxation rates. Vibrational transitions (predominantly changes of one 
vibrational level) are found to be accompanied by simultaneous rotational 
transitions (predominantly changes of zero and two rotational levels) with 
the result that vibrational transition rates are increased by 50 percent or 
more . 


The three-dimensional theory enables one to determine both the gradient 
and the magnitude of the potential, whereas only the gradient can be deter- 
mined with one-dimensional theory. The theory can be fit to data reasonably 
well by appropriate choice of an ’’effective” interaction potential. This 
potential is considerably steeper and of shorter range than potentials appro- 
priate for scattering. This is consistent with the concept that many inter- 
action potentials exist for molecules, just as for atoms. We conclude that 
the steeper inner potentials are primarily responsible for vibrational transi- 
tions, whereas the outer potentials are primarily responsible for scattering. 


INTRODUCTION 


Rates of vibrational transition produced by molecular collisions are 
involved in a number of gasdynamic problems. Vibrational modes of energy 
influence shock-wave structure, cause dispersion of sound, and are responsible 
for sizable variations in the equation of state, for example. The rate at 
which these modes are excited becomes important whenever it is comparable with 
the rate of change in gas properties generated by flow processes, sound waves, 
or other means. 

Vibrational modes of the atmospheric gases are excited so little at 
normal temperatures that they can often be neglected entirely. However, aero- 
dynamicists grew interested in vibrational relaxation with the development of 
supersonic flight; at that time, gas temperatures sufficiently high to excite 



vibrational energy of air were produced in flow about vehicles and in shock 
tubes. Interest continued to grow as higher speed flight generated still 
higher temperatures and dissociation of molecules became an important energy 
sink in the flow, since it was soon recognized that the rates of dissociation 
and recombination were directly coupled to the vibrational transition rates 
(refs. 1-3). Recently, the gas laser has kindled still further interest in 
vibrational excitation and coupled rotational transitions, since the optical 
transition stimulated in some gas lasers occurs between a select pair of 
vibration- rotation levels, and the rate of population or depopulation of these 
levels by collision is one of the factors affecting laser operation. 

So many studies of vibrational excitation have been made that it is 
hardly practical to recall all of them here. However, the principal methods 
of attack on the problem can be summarized by reference to a few selected 
papers. A classic paper by Landau and Teller (ref. 4), in 1936, analyzed the 
one-dimensional collision excitation of harmonic oscillator vibrations in 
connection with the dispersion of sound. Using simple, classical arguments 
about the form of the impulse produced in collision, they deduced that, to a 
first approximation, the vibrational relaxation rate should vary as 
exp [- (0/T) ^ . Probably every set of measured vibrational relaxation rate 
data ever obtained for diatomic gases has been compared with their result, 
generally with reasonably good agreement (see refs. 5-9, for example). The 
Landau-Teller model was refined somewhat by Bethe and Teller (ref. 10) in 
1945. In 1952, Schwartz, Slawsky, and Herzfeld (ref. 11) carried through a 
formal quantum treatment of energy exchange during one-dimensional collision 
between harmonic oscillators, based on the work of Zener (refs. 12 and 13) and 
Jackson and Mott (ref. 14). Takayanagi (ref. 15) independently published 
similar results. This work has been widely accepted as the most rigorous 
practical analytic approach to the problem. Schwartz and Herzfeld followed 
their study with a quantum treatment of three-dimensional collisions (ref. 16) 
based on the method of partial waves (ref. 17) in a manner similar to that of 
Takayanagi (ref. 15) . Although they did not derive expressions for the cross 
sections - which, in principle, can be done with three-dimensional theory - 
they were able to conclude that the transition probabilities have the same 
functional form in the three-dimensional case as for the one-dimensional 
collisions. Treanor, Rich, and Rehm (ref. 18) used the one-dimensional 
Schwartz-Slawsky-Herzfeld theory to investigate effects of anharmonicity on 
vibration-vibration exchange. 

An alternative approach is the semiclassical or impact parameter method 
in which the classical trajectory of motion between gas particles is used to 
obtain a time dependent perturbation potential; then transitions produced by 
this perturbation are calculated by quantum principles. Since the de Broglie 
wavelengths of heaA^^ gas particles are normally much smaller than the scale- 
size of the potentials involved, the classical trajectory is a good approxima- 
tion and the semiclassical results are potentially as accurate as needed for 
most practical purposes. In fact, Rapp (ref. 19) showed that a complete 
classical theory for the one-dimensional collision excitation of harmonic 
oscillators, from the ground state to the first excited level, gives exactly 
the same results as the complete quantum treatment by Herzfeld (ref. 20) . 
However, in the more general case involving transitions between upper vibra- 
tional levels, or vibration-vibration exchange, it is necessary to evaluate 
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the wave function transition matrix elements to get correct results. Rapp and 
Sharp (refs. 21,22) used the semiclassical method to investigate the vibra- 
tional excitation produced by such high energy collisions that small perturba- 
tion theory fails, and Rapp and Golden (ref. 23) analyzed the resonant and 
near resonant vibration exchange with the semiclassical approach, all for one- 
dimensional collisions. Kemer (ref. 24) developed analytic expressions for 
the quantum mechanical transition probabilities of harmonic oscillators, sub- 
ject to large perturbation forces where multiple quantiom jumps occur, as func- 
tions of the energy transfer to a classical forced oscillator; and Treanor 
(refs. 25,26) showed that these solutions are consistent with the numerical 
solutions obtained by Sharp and Rapp (ref. 22) . Of the three basic methods 
available (quantum, semiclassical, classical), the semiclassical method is 
usually simplest, although the classical collision trajectory used in it does 
not correctly conserve total energy in the collision. 

Although one-dimensional collision studies lead to the functional form of 
the transition probabilities, it is necessary to make an educated guess for a 
”steric factor” or a "total cross section” in order to transform such transi- 
tion probabilities into excitation cross sections and ultimately into rate 
coefficients. An absolute derivation of the cross sections must be based on a 
three-dimensional model of the collisions. Calvert and Amme (ref. 27) related 
the one-dimensional treatment to a three-dimensional model in an approximate 
way in order to arrive at a cross section, but did not fully consider the 
three-dimensional details of the collision process with coupled rotations. 

Mies and Shuler (ref. 28) extended Shuler and Zwanzig’s treatment (ref. 29) of 
harmonic oscillator transitions to three dimensions, but for impulsive colli- 
sions that are recognized to be not very realistic, and by numerical methods 
that do not yield analytic relations. Thus, there is need for a further look 
at three-dimensional collision effects on vibrational transition. 

The present paper extends the semiclassical method to three-dimensions 
and provides expressions for the vibrational transition cross sections and 
rate coefficients, including the effects of simultaneous rotational transi- 
tions, in terms of exponential interaction potential parameters. The method 
is then further extended to account for conservation of energy, and the 
effects of distortions in the classical trajectory, caused by reversible and 
irreversible transfer of energy to the vibrational mode, are assessed. The 
values of the potential parameters are found which best fit available experi- 
mental data for vibrational relaxation. Simple analytic approximations are 
deliberately extended as far as possible, even at the sacrifice of some 
accuracy, so that the general functional behavior of the cross sections and 
rate coefficients for vibrational transitions can be followed more readily. 


GENERAL DISCUSSION OF THE PROBLEM 


The problem considered here starts from a given interaction potential 
between colliding molecules and ends with the derivation of a rate coefficient 
and a relaxation time for vibrational transitions. The analysis proceeds by 
the following steps: 
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The interaction potential between a homogeneous diatomic molecule and an 
inert collision partner is assumed to be a simple sxiperposition of exponentials 

U = A E (1) 

i 

where the rj^ are the distances between the inert gas particle and each atom 
of the target molecule. The constant A establishes the energy scale, and 
L is a characteristic length that establishes the potential gradients. In 
the more general case, where the vibrationally excitable molecule is inhomo- 
geneous, separate constants Aj[ and must be used for each atom of the tar- 
get molecule, of course. According to results from scattering and viscosity 
measurements , or from quantum mechanical calculations - for the few cases 
where they are available - it appears that typical values of A are the order 
of several hundred eV, and that L is a fraction of an Angstrom. 


The classical trajectory of relative motion between the gas particles is 
next determined for the spherically symmetric part of the interaction poten- 
tial. This trajectory defines the impulse-like perturbation of the target 
molecule; the perturbation U is a function of time t, the coordinates of 
the molecule q (which represents a set of generalized coordinates q^) , the 
collision energy E, the impact parameter b (initial miss distance between 
the approaching gas particles), as well as the potential parameters A and L. 

U = U(t,q,E,b,A,L) (2) 


If the perturbation is small, or of very short duration, the probability 
of transition from state n to n* may be calculated according to small- 
perturbation-quantum theory. The calculation involves the integration of U, 
weighted by the product of the initial and final state wave functions, over 
all space and time. Thus, the variables t and q are eliminated, and the 
transition probability takes the form 


P 


nn' 


P 

nn 


(E,b,A,L) 


(3) 


Generally, small perturbation theoiy will suffice when comparisons can be made 
with experimental vibrational relaxation times. However, if the analysis must 
be extended to very high collision energies, the small perturbation result 
exceeds unity and the effects of coupling between adjacent states must be 
taken into account. With the approximations used by Kemer (ref. 24), the 
multiple transition probabilities that result from such large perturbations 
can be expressed as functions of the small perturbation transition rate from 
the ground state to the first excited state. Thus, in one form or another, 
the transition probabilities have a functional dependence as expressed in 
equation (3) . 


The collision cross section for transition is just the sum of the 
elements of cross section 27rb db weighted with their corresponding transi- 
tion probabilities. Thus, cross section is a function of E, A, and L only 

Snn,(E,A,L) = JT” db (4) 
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The rate coefficient is obtained by averaging the cross section over all 
collision energies. In most cases, collision energies have a Boltzmann dis- 
tribution, and the rate coefficient may be expressed as 


a^j^,CT,A,L) = u jr“ (x,T,A,L)x e""" dx (5) 

The variable of integration x is the dimensionless collision energy E/kT 
and the factor u is the root mean square velocity. The dependence of cross 
section on E is transformed to the dependence on x and T for the purpose 
of integrating over all x, so the rate coefficient appears as a function of 
temperature T as well as the potential parameters A and L. In the more 
general case, the potential could be a function with as many parameters as 
desired, of course. 

Finally, if the molecules are assumed to be harmonic oscillators, the 
relaxation time is very simply related to the rate coefficient (ref. 20) 

'^(T,A,L) ~ naio[l - exp (-ho)/kT) ] 


where n is the molecule number density. 

With this introduction to the steps to be followed, we are now prepared 
to discuss the semiclassical model of vibrational excitation in more detail. 


SEMICLASSICAL TRANSITION RATE MODEL 


Quantum Perturbation Theory 

A brief review of the quantum mechanical perturbation method for 
calculating transition rates will be helpful for the purpose of discussing the 
perturbation potentials. The method begins by expressing the wave function 
'JP of the system during collision as a linear supeirposition of the orthonormal 
set of wave functions that describe the system in its unperturbed states 

n 

The generalized coordinate q is londerstood to represent all the space 
coordinates involved in the problem, and is the energy of the nth 

quantum state of the unperturbed system. If the total wave function is nor- 
malized, the squares of the absolute values of the coefficients anCt) can be 
interpreted as the probability that the system will be in the state n at the 
time t. Thus, these squares sum to unity. 
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1 


( 8 ) 


2 

n 

The wave fimction of equation (7) is introduced into Schroedinger ’s time 
dependent equation 


Hy = ih (9) 

in which the Hamiltonian operator H is expressed as the sum of the 
Hamiltonian for the unperturbed system H^ (which is independent of time) and 
the perturbation energy U created by collision (which is a function of both 
molecular coordinates and time) 

H = H°(q) + U(q,t) (10) 

When both sides of equation (9) are multiplied by the complex conjugate wave 
functions exp(iEnit/h) and integrated over all space, the differential 

equations for the coefficients an» are obtained 

^n' ^ ~ (t)e'-“nn't ^ . . . (11) 

^ n 

where the circular frequency is 


nn ' 


En' - En 
h 


( 12 ) 


and the transition matrix element H , is 

nn’ 




(13) 


Equations (11) are a coupled set that must be solved simultaneously for 
all the coefficients a^, subject to the given initial conditions at t = -00. 
In the present case, the molecule is known to be in the state n before 
collision, so the initial conditions are 


( 1 if n' = n 

0 if n' 7^ n 


(14) 


I 1 2 

The final values |a ,(+«») | represent the probability that the collision has 
caused transition to the state n*. 
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In some cases the deviation from the initial state is small throughout 
the collision. Then the differential equations for the coefficients a , are 
uncoupled and become 




a„, = - i H , 
ri' nn' ^ 


(15) 


The transition rate P , is 

nn’ 


— nri 


(16) 


The quantity ^nn* interpreted as a transition probability only 

if it is small compared to unity. This will normally be the case at the rela- 
tively low collision energies which obtain when kT << In fact, over 

the range of temperatures where measured relaxation times are available for 
comparison with theory, P^^^? is usually less than unity, even where kT is 
somewhat larger than hoj. However, at very high temperatures, P^in’ 
exceed unity and cannot then be identified as a probability. Physically, it 
represents a transition rate between states in the coupled set of equa- 
tions (11). In principle, the coupling can be taken into account for strongly 
perturbed harmonic oscillator transitions, by an approximate method developed 
by Kemer (ref. 24) , which yields transition probabilities as a function of 
Pqi. Thus, calculation of Pqi is the first step in the solution of the 
problem in either case - weak or strong coupling. 

The matrix element needed for the calculation of Pnn» 

obtained from equation (13) once the perturbation energy U(q,t) is defined. 

An approximation to this perturbation energy for homonuclear, diatomic 
molecules is developed in the next section. 


Approximate Perturbation Energy for Homonuclear 
Diatomic Molecules 

The coordinate notation used to analyze the three-dimensional collision 
between a homonuclear, diatomic molecule and an inert collision partner is 
shown in figure 1. The center of the diatomic molecule is the origin (only 
one atom of the molecule is shown in figure 1; the second atom is understood 
to be symmetrically placed about the origin). The trajectory of relative 
motion between the particles is assumed to be determined by the spherically 
symmetric part of the interaction potential, in which case the trajectory 
(shown by the dashed curve) lies in a single plane (the x-y plane of fig. 1) . 
The point of closest approach is taken to be on the x-axis a distance a 
from the origin. The distances between the inert particle and the center of 
the molecule, the center of atom 1, and the center of atom 2 are designated by 
r, ri, and r 2 , respectively. The angle between the x-axis and the vector r 
is X. The distance between atoms 1 and 2 is p, and the orientation of the 
molecular axis is defined by the polar angle 0 from the z-axis and the 
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Figure 1.- Diagram showing collision coordinates. 


azimuthal angle ^ from the x-axis. Then the distances ri and X 2 are 

ri2- = r^-jl - ^ sin 6 cosC<l> - X) + ] 

2 

X 2 ^ = r^ Fl + — sin 0 cosC4> “ X) -t 1 



Collisions are assumed to be weak enough that (p/r) < 1 for all r. Then, 
reasonable approximations for the distances ri and T 2 are 

ri - r - y sin 0 cos(<j> - X) (18a) 

r 2 - r + Y sin 0 cos(<l) - X) (18b) 

At this point the perturbation potential is taken to be a linear super- 
position of exponential repulsions, as in equation (1) 


U = - 2Ae"^^^ coshj^^ sin 0 cos (cj) - X)j (19) 

The bond length p is very nearly equal to its equilibrium value p^, 
and the angle X is small in the region where the perturbation contributes 
appreciably to the Fourier transform. Accordingly, the perturbation is 
expanded to terms of first order in p - Qq and x 


D e 

sin 6 cosC<f)- X} - — sin 


2L 


cos 


Pe . 
^ 2L 


0 sin 


2L 


sin 0 cos 


( 20 ) 


Then, if 6 = Pg/2L 

U - 2Ae"^^^ jcosh (6 sin 0 cos ((>) + X6 sin 0 sin c}) sinh(6 sin 0 cos cf)) 


P-Pe 

2L 


sin 0 cos 4> sinh(6 sin 0 cos <j))J 


( 21 ) 


The spherically symmetric potential used to determine the collision 
trajectory is obtained by averaging over all configurations. The last two 
terms vanish because of asymmetry in (p and (p - p^) , and the part remaining 
is 


U - 2Ae’^^^ cosh (6 sin 0 cos cj)) = 2Ae 


( 22 ) 


This spherical averaging procedure is common practice in quantum treatments, 
and often seems justified by the results. It is also justified in classical 
treatments of the motion when a >> p. Although a may not exceed p by a 
large margin, the inequality a > p is usually satisfied at least, so the 
error in the approximation of spherical symmetry does not grow unreasonably 
large. 
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The constant A will be defined as 


• A 


2A 


sinh 6 
6 


(23) 


This constant represents the effective potential size parameter for the total 
spherically symmetric potential when the effects of the two separate atoms in 
the molecule are added. The perturbation is now expressed as 




sinh 6 


cosh (6 sin 0 cos cf)) + X6 


sin 0 sin ^ sinh(6 sin 0 cos (0 


+ 


P-Pe 

2L 


sin 0 cos (() sinh (6 sin 0 cos 



C24) 


where = Ae the spherically averaged potential at the point of closest 

approach . 

The first two terms in equation (24) are functions of only the rotational 
coordinates and are responsible for elastic scattering (aZ = 0) and pure 
rotational transitions (AZ ^ 0) . The third term is linear in (p - p^) and 
leads to vibrational transitions as well, which in the harmonic-oscillator, 
small '♦perturb at ion approximation consist of single vibrational quantum jumps 
(Av = ±1). However, the appearance of the factor sin 0 cos cf) sinh(6 sin 0 cos (\>) 
in this term means that these vibrational transitions can be accompanied by 
simultaneous rotational transitions, a feature missing from the one- dimensional 
collision model. 


As far as vibrational transitions are concerned, the effective perturba- 
tion is just the third term, which may be expressed as 

CO 

II = Upg" / P'Pe\ \^ (6 sin 6 cos . . 

sinh 5 V 2L /X (2n - 1)! 

n=i 

From the symmetry of the angular coordinates, rotational transitions are 
limited to an even number of quantum jumps. Thus, the symmetry of the 
diatomic molecule’s rotational wave function is preserved in these transitions, 
consi'Stent with the expectation that collision perturbations are not likely to 
change nuclear spin. If nuclear spin is unchanged, rotational symmetry must 
be preserved to maintain antisymmetry of the total wave function. 

The summation of equation (25) is uniformly convergent, even though 6 
is normally greater than unity. For typical values Qq = I ^ and L = 0.2 A, 

6 is about 2.5. Then the first two terms of the perturbation are about 
equal in mangitude, but the third term is only 15 percent as large as the 
first two, while the fourth term is only 2 percent as large and succeeding 
terms become very small. In view of these magnitudes, it is obviously neces- 
sary to carry at least the first two terms of the expansion; to simplify the 
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analysis we shall truncate the perturbation at this point, however. 




62 

sinh 6 



6 cos 2 



sin^ 


0 cos'* 


♦) 


( 26 ) 


In this case, the coupled rotational transitions are limited to AZ=0, ±2, ±4 
The various overlap integrals and the matrix elements resulting from the per- 
turbation of equation (26) are worked out in appendix A. The truncation does 
not seriously affect the value of sinh 6, so this factor is retained without 
change. 

The factor exp [- (r-a) /L] contains the functional dependence on time 
which determines the magnitude of the Fourier transforms in equation (16). 

The classical equations of motion will next be used to develop an 
approximation for this time -dependent factor. 


Classical Trajectory Impulse Functions 


The classical trajectory for a spherically symmetric potential lies in a 


single plane as shown in figure 2 (the 


X 



Figure 2.- Classical trajectory of relative motion 
between particle P and molecule M, 


x-y plane of fig. 1) . The velocity 
of relative motion at large r is u, 
and the impact parameter b is the 
initial miss distance. Energy and 
angular momentum changes involved in 
the internal state transitions of the 
diatomic molecule will be neglected 
compared with the kinetic energy and 
angular momentum of the collision 
partners, so that the latter can be 
treated as conserved quantities. 

(Later some adjustments will be made 
to account for conservation of total 
energy.) Conservation of the angular 
momentum of the collision partners 
requires that 


mr^X = mub 


(27) 


while conservation of energy requires that 


^ (f^ + r^X^) + U(r) 


mu^ 

2 


(28) 


where m is the reduced mass of the two collision partners and U(r) is the 
spherically symmetric average potential obtained in equation (22) . Solving 
for X from equation (27) and substituting in equation (28) , one obtains 
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( 29 ) 



2U(r) _ 
mu^ r^J 


At the distance o£ closest approach r vanishes by symmetry (for an elastic 
collision) and from equation (29) 





(30) 


Precise knowledge of the potential is most important near the turning 
point where the perturbation is largest. Expansions about the point of 
closest approach (t = 0) yield the even function for the perturbation impulse 

TJ(t) 


u(t) 


y -[(r-a)/L] 


^ y^g-[Crot^/2L) + (r^^t'+/24L) + . . .] 


(31) 


and the odd function for the angle X(t) 


1 3 

X(t) = X^t + Xq . . . (32) 

All time derivatives at the turning point may be evaluated in terms of L, b, 
and a (i . e. , rQ) 


to = 


^0 u^b^ 
mL ^ rt3 


„iv 


/ ^o 3u^b^\ •• 


C33) 


x„ = 


ub 




6 * ^)' 

\ mu^Lb^/^ 


(34) 


For most purposes it would be sufficient to use just the first term in 
the expansion for X, but, as we have seen, this quantity is not needed for 
calculations of vibrational transitions anyway. Only the exponential function 
in equation (31) is needed for this purpose, and this function is very sensi- 
tive to the exact value of r. The first exponential factor falls too rapidly 
as t^ increases. The following factors oscillate in sign of the exponent, 
and the total expression does not converge very rapidly. Thus equation (31) 


12 



is not very useful as it stands. However, for the exponential potential form 
an exact^solution does exist for head-on collisions (refs. 12, 13, 19) 

(b = 0, Uq = mu^/2) 


U(t) = Uq sech^Cat) 


(35) 


The constant a is u/2L when b = 0. Now it is anticipated that the most 
effective collisions are those where the miss distance is small, since these 
penetrate deepest into the potential and cause the strongest perturbations. 
Therefore, we choose a solution with the same form as equation (35) that will 
reduce to the exact solution when b = 0, but let the constant a take the 
value that gives the correct first-order terms in the expansion of the three- 
dimensional trajectory about the point of closest approach 


2 _ _ 1 /^o u^b^\ 

- 2L “ 2L \mL J 


(36) 


We shall define e as the fraction of the collision kinetic energy that is 
transformed into potential energy at the point of closest approach 


£ 



(37) 


In terms of e, the constant a may be expressed 





(38) 


This value of a gives the correct form near the turning point at the expense 
of some mismatch in the asymptotic behavior. The particles should recede from 
one another with a velocity u at very large time, while the asymptotic form 
of equation (35) 


U(t) 


... — 2 at 


4e 



-(ut/L){l-(l-£)[l-(2L/a)]}^'^ 

e 


(39) 


1/2 

corresponds to an escape velocity u{l - (1 - c) [1 - (2L/a)]} . Even so, 

equation (35) is a far better approximation than the simple expansion of 
equation (31), for the asymptotic behavior is correctly described as £ ^ 1. 
Actually, we need not be particularly concerned with the slowly varying per- 
turbation far from the turning point, since this portion does not contribute 
appreciably to the Fourier transform in equation (16) . Rather, an approxima- 
tion is required which is nearly correct in the region where the rate of 
change dU/dt has passed its maximum, as this will pick up the major contribu- 
tion to the Fourier transform. Equation (35) should be a reasonable approxi- 
mation in this sense as long as e - 1, and this may be sufficient to pick up 
the important contribution to the vibrational transition cross section. 
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The perturbation potential of equation (26) is next used to calculate the 
transition probabilities, the cross sections, and the rate coefficients for 
vibrational transitions produced by low energy collisions, where the small 
perturbation method is valid. These calculations will illustrate the proce- 
dures outlined in the General Discussion of the Problem. 


TRANSITION PROBABILITIES, CROSS SECTIONS, AND RATE COEFFICIENTS 


Vibrational Transition Probability 

That part of the perturbation fxmction which contributes to vibrational 
transition has been approximated 

u = e sech^ (at) ® ® '*’) 

The perturbation matrix element for simultaneous transition from quantum 
state (v, Z, m) to state (v', Z', m') is 

H(v,Z,m;v',Z',m') = B(v,Z,m;v*,Z' ,m') sech^ (at) (41) 


where the constant 3(v, Z, m; v', Z', m') is the overlap integral 
3(v, Z,m;v' ,Z> ,m’) = ® 4- +^sin^ 6 cos'* 


vZ-m/ 


(42) 


These constants are evaluated in appendix A for the usual rigid-rotator, 
harmonic-oscillator wave functions. Changes in vibrational quantum number are 
limited to Av = ±1, while changes in rotational quantum number are limited to 
b.1 = 0, ±2, ±4. The changes in magnetic quantum number are limited to Am=0, 
±2, ±4 also, and these correspond to reorientation of the angular momentum 
vector in space as result of the collision. The energy of the rigid rotator 
does not depend on magnetic quantum number, so the transitions involving dif- 
ferent initial values of m can be averaged to obtain a total transition 
probability for (v, V) (v', Z’). Since the probability is proportional to 
the matrix element squared, the values of 3^ must be averaged to get the 
total matrix element squared for a given transition (v, V) to (v+Av, l+Al) . 

3^(Av,AZ)= X j ^ 3^(v,Z,m;v+Av,Z+AZ,m+Am) (43) 

Am=0,±2,±4 ^ m J 

All of the transitions involved in this average have the same value of 
circular frequency 
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0 ) = (jd(Av,AZ.) = 


h 


(44) 




Although it is difficult to average over m in general, the average 
approaches a simple limiting value for t » 1, as shown in appendix A. This 
limit will be used here, inasmuch as t is the order of 10 or greater for 
most of the collisions that occur in gases at the temperatures of interest, 
kT >=« hw. Furthermore, the limiting average for large t is independent of 
t, which is fortunate in view of the averaging over a Boltzmann distribution 
of rotational states that must be performed eventually. 

According to the results in appendix A, the final values for the averaged 
matrix elements squared are 


B^C+1,0) 


b2 (±1,±2) 

B2(±i,±4) 


^ \ r v-H (1/2) ± (1/2) 1 

Vl5 35 630Asinh2 6/L SL^ym/h J 


/I , 52 ^ 26V\| 

f \ 

v+(l/2)±(l/2)' 

V30 105 2835^' 

(sinh2 

8L^pco/h 


( ^ 

/ 6“* \ 

rv+(l/2)±(l/2)~ 

^22680 y 

^sinh^ sj 

8L2yo)/h 


(45a) 


(45b) 


(45c) 


where p is the reduced mass of the diatomic molecule. The next term in the 
expansion of the perturbation function of equation (25) will contribute some- 
thing to the terms of order 6^, but not to the lower order terms. 

Note that diatomic molecules with electronic states other than should, 
strictly speaking, be described by symmetric top wave functions. However, the 
projection of electronic spin and angular momenta on the internuclear axis is 
rarely more than 3/2 quantum units in molecules of interest, so at high rota- 
tional quantum numbers Z, the perturbation of rigid- rotator wave functions by 
electronic momentum should be small. Accordingly, we apply the above rigid- 
rotator matrix elements to diatomic molecules in general. 


The Fourier transform of the sech^(at) is a well-known integral 


£ 


(sech^ at)e^^^ dt = 


TTO) 

72 


cs ch 


TT(jl) 


27TOJ -(77ca/2a) 


2a 


(a)/a) >>1 


(46) 


Thus, in the limit as (w/a) >> 1, which is a condition for small perturbation 
theory to be valid, the transition probability becomes 


"»2 

o c /mu^\ 27UJL) -(TTco/a) 

® h \~r-) ir] 


(47) 
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where the indices indicating the particular transition involved are dropped 
for economy in notation. The terms P, m, and 6 are understood to represent 
the values of P(Av,Al), (o(Av,AZ), and 3 (Av,AZ). The latter two constants are 
given by equations (44) and (45) , respectively . 


Vibrational Transition Cross Section 
The collision cross section for the transition may be expressed 


S = Po f 2Trb db (48) 

J Q O 

where Pq represents the transition rate for head-on collisions when b = 0. 
According to equation (47) 



/aQ\^ " t C^o/^) ~ 


and from equation (39) 



a 


1 


(1 



(49) 


(50) 


where a^ = u/2L, o is the distance of closest approach, and will be 

used to designate that particular value of a which occurs for head-on colli- 
sions, that is, when b = 0. From the relation between b and a (see 
eq. (30)) 

b^ = a2(l - e) (51) 


it follows that 


2.b db . -™2[l - 


For the exponential form of the interaction potential assumed here 

-[(o-Oo)/L] L 

^ ® ’ de e 


and in this case, equation (52) becomes 


2irb db = -tto^^I + 


(52) 


(53) 


(54) 
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Substituting equations (50) and (54) into the pre-exponential factor of the 
integrand in equation (48) , one obtains 


S 



1 


TTO 


2 


+ 


a 



- (TToi/ao) [Cao/a)-l] 
e 


de 


It can be anticipated that if TTw/a^ >> 1, the major contribution 
integral will come from the region where e - 1 and a - a^. Therefore 
integrand of equation (55) will be expanded to terms of first order in 
From equation (53) 


(!-£)+...] 

Thus 


TTa' 

1 

2 

TTO ^ 
O 

1 + (1- e) + . . . 

L °o J 

1-^1 

o 

'' 1-e'^ 

^ e J 

' 

1 

2L ^ 

■ — (1 - o + . . . 

°o J 



while from equation (50) 


= c^o(l - k 1 




Thus, to this approximation, the cross section becomes 


S = TrOo^p^ 


I 


1 - (iTw/ 2 ao) [1 - (^L/ct ) ] (1-e) 


de 


This integration is easily performed to give 


^^0^0 j C^L/Oq) ] I 

' " fe)(^ - a ' ' 

In the limit as (7ra3/2ao) >> 1, the cross section takes the form 


TTOo^ 

2 

TTO^ 

P — ^ 

/23mu^\ /ttoj \ -(Troo/an) 

- S) 

o - OT 

1 1 - — 
’ Oo 

( hao ) baj ® 


(55) 

to this 
the 

( 1 - 0 . 

(56) 

(57) 

(58) 

(59) 

(60) 
(61) 
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The factor (7ro)/2ao) [1 - (2L/0Q) ] may be thought of as the effective 

total cross section for vibrational transition to be used with the one- 
dimensional collision transition probability Pq. It can be seen that this 
steric factor is not constant, as usually assumed, but because of the three- 
dimensional nature of the collision, is approximately proportional to the 
velocity, u, as suggested by Calvert and Amme (ref. 27). 

In the other limit as 7ro)/2aQ[l - (2L/aQ)] << 1, the effective cross 
section approaches according to equation (60). This is a physically 

plausible result, though the assumptions made tend to break down in this limit. 


Vibrational Transition Rate Coefficients 

The cross section is next averaged over a Boltzmann distribution of 
collision energies to get the rate coefficient a 

a = ^ Sxe ^ dx (62) 


where x is the dimensionless collision energy mu^/ (2kT) , m is the reduced 
mass of the collision partners, and u is the mean velocity (8kT/7rm) ^ . The 

symmetry number (s) is unity for collisions between unlike molecules, but must 
be 2 for collisions between like molecules, in order to avoid counting 
collision systems twice in this case. 

Some allowance must now be made to account for conservation of energy in 
those collisions that cause vibrational transition. Rapp (ref. 19) lets the 
effective collision velocity u be the average of initial and final veloc- 
ities (ui + Uf)/2, and Herzfeld (ref. 20) points out that this substitution is 
necessary to reconcile the classical and quantum results. To terms of second 
order this is equivalent to letting the effective collision energy be the 
average of the initial and final energies. 


mu2 _ ^ hco 

2 2 “ 2 


(63) 


where is the initial velocity. The use of the ± sign depends on 

whether the transition is to the adjacent lower or upper vibrational state, 
respectively. The effective dimensionless collision energy x is, 
accordingly 


X = Xi 


^ h(jL) 

' 2idf 


(64) 


where Xj^ is mui^/2kT, the initial value of x. The cross section has been 
derived as a function of x, while the Boltzmann distribution of collision 
energies given in equation (62) is really a function of x^, the initial 
value. The variable of integration is changed from Xj[ to x, and the rate 
coefficient becomes 
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a 


ij ± fiitJL)/ 2kT) 
— e 
s 


/" ^ 

*'hu/2kT 


(x)^x + 


hco \ 

2kT) 


“X 


e 


dx 


(65) 


The factor exp(±ho)/2kT) is exactly that required so that the ratio of the 
forward rate to the reverse rate is the Arrhenius factor exp(-ha)/kT) , which 
preserves detailed balancing at equilibrium. The lower limit of the integral 
is found to be the same for either excitation or deexcitation, since the exci- 
tation cross section vanishes for collision energies less than ho). 


The cross section has been found as a function of the effective velocity 
u in equation (61) 


S 



16 TTg^m^tjL)L^u 


- (2tTCl)L/u) 
e 


(66) 


It will be convenient to express this result in terms of the dimensionless 
collision energy x and a characteristic energy 


Ec 


moj^L^ 

2~ 


(67) 


Then equation (66) becomes 


S 



( 68 ) 


where the constants F and G are 


/Ec3kT\l/2 
F - 64,6^ ( ^ 


(69) 


, , 1/2 

■ - (w) 


(70) 


Substituting equation (68) into (65) , one obtains for the rate coefficient 


±(hw/2kT) 

a = — e 
s 


r 1 / 2 / \ 

^ ^ If ■ 2kT/ 

hu/2kT ^ 


-[x-h(G/x1/2)] 

e dx (71) 
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The integrand has a sharp maximum primarily determined by the factor 
exp- [x - (G/x^/^)]. This factor has its maximum where the exponent is a 
maximum. 


^m = 




( 72 ) 


The characteristic energy Ej, is several eV for most diatomic molecules, 
and the constant G is the order of 10 to 100; therefore, the maximum occurs 
at rather large values of x, just as we assumed when the corrections to the 
pre-exponential factors and to the lower limit of the rate coefficient integral 
were ignored. 

The pre-exponential factors x^^^[x - (ho)/2kT)] and in equation (71) 
are evaluated at x^, and the exponent is expanded about x^ 


1/2/ _ 

^m (^^m 2kT / 


3/2 


m 


(73) 




(74) 


TTOf 


a = — 
s 


a - 


2L ^m 


-3x„,± (ho)/2kT) 


U| "-0 1 p^J/2 Q f 


“(3/^Xni)(x-Xin)^ , 

e dx 


Or 


h(D/2kT 


m 


The integral in equation (75) is approximately 


/: 


-3y2/>^Xm 
e dy 


/4ttx 


m 


(75) 


(76) 


and the constant F expressed in terms of Xj^ is 

2/8kT\ 9/2 


F = e 

so the rate coefficient finally becomes 


/8kT\ 

Ivho) j ^ 


(77) 


a = 


s 


l/2/_ 2 


/4tt\ /^®o \ ^/8kT\ „ 1 

( 3 ) L 2L ) ^ 


3/2 g 3x^-(*^<^/2kT) 


m 


(78) 
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where the distance of closest approach for head-on collisions at the energy 
Xj^kT is 


Oq = -L Zn 


^mui 


TO 


2A 


= L Zn 


(v^t) 


(79) 


The rate coefficient, expressed in terms of temperature and the 
characteristic energy Ec> is 


3' 


\/8kT - sCtt^Ec/IcT) ^ '^±ho)/2kT 

vrj L ^ l\Trm j V ho) / Wt / ® 


(80) 


TO 


where the total cross section at the velocity Uj^ is 


fra. 


ttL2 

^ TT^k^f^EcV 


o. 


(81) 


1 - 


m Zn(A3/7r2k2T2Ec) 


Equation (80) demonstrates the familiar result (refs. 4, 11, 19) that Zn a 
varies essentially as 

Rapp (ref. 19) called attention to one additional factor, due to long- 
range attractive forces in the intermolecular potential, which can become 
important at low temperature. If the depth of the potential well is A, the 
collision energy is effectively increased by A before it is absorbed in 
climbing the potential barrier. This amounts to a shift in the Boltzmann 
distribution of collision energies by the factor exp(A/kT). Since this fac- 
tor does not enter the integrations, it simply multiplies the final results of 
equations (78) and (80) , 

Equation (78) or (80) gives the rate coefficient for a specific vibration- 
rotation transition. The values of and 03 will be different for different 
transitions; recall that the indices specifying the exact transition were 
omitted for convenience. Now, however, the total rate coefficient for a 
given vibrational transition, a(Av), is obtained by summing the coefficients 
for all possible rotational transitions, a(Av,AZ), and averaging over the 
initial distribution of rotational states 


a(Av) = a(Av,AZ) (82) 

AZ=0,±2,±4 

This is accomplished by substituting for 3^ in equation (80) an appropriate 
weighted average, which is derived in the next section. 
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Effect of Rotational Excitation 


To a first approximation, the energies and the circular frequencies 
oj(Av,AZ) are independent of the rotational quantum number Z. In the limit as 
Z >> 1, the matrix elements averaged over all values of m, 3 ^(Av,AZ), are 
also independent of Z, as shown in appendix A. Thus, the final matrix ele- 
ment for a given vibrational transition Av may simply be taken as the sum 


e^CAvD 


E 

AZ=0 ,±2,±4 


32(Av,AZ3 


(83) 


and this value used in equation (80) will give the first-order approximation 
to the total rate coefficient for a given vibrational transition v to 
V + Av. 

In the second approximation, the small variations of oj are taken into 
account. To first order in v and Z, the energy of a given state is 


E = ho)(v + + 1) 


(84) 


where B is the rotational energy constant. The circular frequencies for the 
allowed transitions are 

,,, ^v+t\v,l+Ll " ^v,l r. 2BAZ /, AZ+1M 
a)(Av,AZ) = i ^ = cjIAv + (3 + -J~]\ (85) 

Since I is normally large compared with unity, [(AZ-+1)/2] will be neglected 
compared with I, and the approximate expression lased for the circular 
frequency is 


)(Av,AZ-) to Av + 


( 86 ) 


Then, approximate corrections to the first-order rate coefficient a(Av,0) are 


f. .7^ r /-A -3Xn,| [to(Av,AZ-)/a3(Av,0) 1^^ ^-1} 

Using equation (86), and setting to(Av,0) = coq, it is foimd that 


( 87 ) 


gj(Av,AZ-) _ 2 7. / 2 BAZ- \ ^ 
Wo ~ \fiwo / 
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Since (26/1103^) ?=« 2x10"^ for typical molecules, the exponent in equation (87) 
may be expanded and 


R(Av,A7) 


"l + / 2BAZ \ ^-Xm(4BAi/ha)o)^ 


(88) 


Dependence of R(Av,A7) on I is removed when the correction factor is 
averaged over a Boltzmann distribution of initial rotational quantum numbers 


R(Av,AZ) 


j|r^^R(Av,AZ) C2Z+l)e 
1=0 


(B/kT) 1(1+1) 


(89) 


An analytic approximation for this average can be obtained by replacing the 
summation with an integral. Set 

, _ 2B|a^| 

ho)o V B j 

g = 


Then 


y2 =^Ut ^ 1) 

7 /1cTV^2 

>'(r) 


R(Av,AZ.) 



(1 + xyy^'^ 


-Cy+g)^ 

e 


2y 


dy 


R(Av, AZ.) 



(1 - Ay) 7/ 3 


-( y - g )^ 


2y 


dy 


\AV 


\Av 


“) 

“) 


(90a) 


(90b) 


The integrand in equation (90a) has a maximum near y = 0 and in equa- 
tion (90b) the integrand has a maximum near y = g. Using these values for y 
in the pre -exponential term, and integrating the resulting simplified integral 
gives 

R(Av,AZ) = 1 - /iT g e® (1 - erf[g]) > 0^ (91a) 

R(Av,Al) = (1 - >^g)^'^^|^l+ /iT g e® (1 + erf[g])j < oj (91b) 
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R(AV,lAll) R(AV,IAl!J 


These are combined to give a total correction factor for a given |aI|. 


R(Av, |aI]) = 1 - /it g e® (1 - erf [g])+ fl - Ag) ^ + /ir g e^ (1 + erf [g])J 


( 92 ) 


5.5- 
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Figure 3.- Correction factor R(Av,|a 3|) for 
I All = 2 and L= 0.2 and L= 0.25. 

?4 — 

N2 L = .25 a 


The correction factors given by 
equation (92) are shown as functions 
of tenperature in figures 3 and 4 for 
diatomic molecules N 2 and 0^ and for 
typical values of L = 0.2 and 0.25 A. 
It is seen that transitions v->v ±1, 
1-^+2 can increase in probability by 
factors of about 4 and that transi- 
tions V V ±1, 1 1 +A can increase 

by factors of 10. The matrix elements 
depend only on the absolute magni- 
tude of the change in Z, so the total 
correction for a given change [AZj 
is the sum of equations (91a) and 
(91b). For all transitions involving 
a given change in vibrational quantum 
number Av, the weighted average value 
to be used for 6^ in calculating the 
total rate coefficient is 


22-/ 


20- 
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N 
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L = .25 A 


\ 



B^(Av) = E R(Av,AI)b2(Av,AZ) 

AZ=0,±2^±4 


(93) 


and 


^ ^ R(Av,AZ-) (AVjAZ.) 
g AZ-=0 ,±2 ,±4 

AZ=0,±2,±4 


(94) 


4*- 

2000 


Figure 


This ratio is shown for O 2 and N 2 
in figure 5 for L = 0.20 and 0.25 A. 
The total transition rate is seen to 
, , I I I I increase by factors of 1.5 to 2 when 

3000 4000 5000 6000 7000 8000 coupled rotations are taken into 

^ account. Such an effect could simply 

4.- Correction factor ^avJa^IT for be swallowed by the empirical cross- 
|az| = 4 and L=o .2 and L=o. 25 . Section factor in the one-dimensional 
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2.1 ■ 


collision model, of course, but in the 
absolute three-dimensional model a 
factor of this size makes some differ- 
ence in the fit of theory to data. 

It may be helpful to summarize 
here the essential formulas for the 
particular rate coefficient olio> 
which is of primary interest. The 
potential parameters A and L, the 
reduced mass m, and the symmetry num- 
ber s are all determined by the two 
collision partners involved; the cir- 
cular frequency coq, the equilibrium 
interatomic distance p^, the molecu- 
lar reduced mass p, and the rota- 
tional energy constant B are all 
characteristic of the diatomic molecule which is vibrationally de-excited. 

In terms of these parameters and the thermal energy kT, the rate coefficient 
for de-excitation from the first vibrational state is 



Figure 5.- Total rotational correction for 
the rate coefficient, R = 


^10 

s \3 / / \Trm / \ / 


-sCtt^Ec/IcT) ^/3+(hwo/2kT) 


(95) 


where the cross section S is 


. ^ uWeJ 

^ " A 

1 -■ 


In 

and the characteristic energy E^. is 

I, 

Ec = 


(96) 


(97) 


The constant 3fo given by 


_ m hojo 5 ^ 


10 


16p sinh^ 6 


So) ^ iSs) 


( 98 ) 


where 
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- - ■ 2 — 2 

R(-l, |aZ|) = 1 - /iT g e® (1 - er£[g])+ (1 - Xg)'^''3|’i + /IT g e® (1 + erf[g])j 


( 99 ) 


with X = (2B|AZ|/ho3o) (kT/B) 1^2 ^nd g = XCir^Ec/kT) 3 ^ 

These results will later be used to make some comparisons with experiment, 
and to deduce some values of the potential parameters A and L which best fit 
vibrational relaxation data. 


CONSERVATION OF TOTAL ENERGY 


We now consider the effects produced when conservation of energy is 
extended to include both the inelastic transfer of energy that occurs when 
vibrational transition occurs, and the elastic but nonadiabatic transfer of 
energy that occurs because the vibrational mode is transiently excited by 
collision (refs. 21, 22). In the following discussion an elastic collision is 
defined as one where no net change of kinetic energy is produced. The elastic 
collision may be adiabatic , meaning that internal modes of energy do not par- 
ticipate in the collision process; or the elastic collision may be nonadia- 
hatio , in which case the internal modes of energy participate transiently. 

The inelastic collision is one where permanent transition occurs in one or 
more modes of internal energy, with a corresponding net change in kinetic 
energy. All inelastic collisions are also nonadiabatic, of course. 

The usual perturbation used in semiclassical and classical analysis of 
the vibrational excitation problem is elastic and adiabatic in form, but, as 
pointed out earlier, the inelastic quality of the collision is approximately 
accounted for by letting the effective collision velocity be the average of 
initial and final velocities. This simple strategem does not fully account 
for inelastic and nonadiabatic effects, however. The perturbation is also 
asymmetrically distorted by the transition, and, in addition, the vibrational 
modes of the diatomic molecule are transiently excited during collision, 
whether or not transition occurs (refs. 21, 22). The corrections required 
because of these inelastic and nonadiabatic collision effects will next be 
assessed. 


Effects of Asymmetrical Distortions to the Collision Trajectory 

When the collision is not adiabatic, the trajectory is asymmetrically 
distorted by permanent transfer of energy from or to the kinetic mode, and the 
classical conservation equation for energy 

E = T + U (100) 

does not hold. In this case, U is a skewed function and such an asymmetric 
function may be represented by 
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( 101 ) 


-(¥) . 


U = Uoe ' “ ' = eEe*^^ sech^ at 

where q is a constant. With this form for the potential, the relative 
velocity between collision partners is 

•^ = 2a tanh at - q 

and the constants a and q may be made to fit the limiting conditions 
r(_oo) = -ui and f(“) = U£. 


(102) 


a = 


Ui + Uf 
4L 

Ui - Uf 
2L 


that 


For the transitions of interest, one quanttun of energy is involved so 


j (ui^ - = ±htjj 


and this equation can be combined with the two previous equations to yield 

,2 _ 


2(ui^ - ^ hoj 

a " (ui + Uf)^ “ 2E 


a = 


(103) 


The probability of a transition is appreciable only when the collision energy, 
E, is large compared with ho), and so, only when q/a << 1. 

With the interaction potential of equation (101), the transition 
probability is 


P* = b" |F*|" 

v,v±l v,v±i \nIo/ I I 


(104) 


where 


F* = 03 


/: 


3^^ sech^ at e^^^ dt 


(105) 


This transform is evaluated in appendix B, where it is found that 

-Tux}/a 


|F*|^ = 4.2(^) (l . ^)e‘ 


(106) 
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Therefore , 




-TTU/a 


(107) 


Except for the factor [1 + (q^/w^)j, this expression for the transition 
probability is the same as the earlier (eq. (47)) result which follows from an 
assumed adiabatic trajectory. This factor (q/w)^ represents the correction 
for skewness in the shape of the interaction potential, and its size is easily 
estimated with the aid of equations (97) and (103) . 


(if - (if (if 


16 EEc 


For collisions of interest, (hto/Ec) < 10"^ and (hw/E) << 1, and so 


(108) 



10-3 


This size for the correction factor leads to the conclusion that the major 
effect of distortions to the trajectory induced by vibrational transitions is 
accounted for by defining (refs. 19, 20) the effective collision velocity, u, 
to be 


u 


u. 


+ 


Uf 


2 


and that the transition probability is only slightly affected by the 
asymmetrical part of the distortion. 


Effects of Adiabatic Symmetrical Distortions 


In order to examine the effect of trajectory distortion caused by 
temporarily stored vibrational energy, it is first necessary to develop an 
expression for this quantity, V, as a fimction of time. This is done in 
appendix C with the help of simple relations, found by Kerner (ref. 24), for 
the energy excited in a classical oscillator and quantum transition probabil- 
ities for a harmonic oscillator. For transitions from the ground state (on 
which interest is here centered for comparison with observed relaxation data) , 
the total number of vibrational quanta excited in all states is given by 
exactly the same expression as by small perturbation quantum theory. 


v(t) = hwP(t) = ^ 



U(t)e^“'*^ 


2 


dt 


(109) 
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In this equation, the subscripts 0, 1 on P and 6 are omitted for convenience, 
and it is assumed that U is an exponential interaction potential. Also in 
appendix C, the coefficients of the expansion of V(t) about the turning point 
(t = 0) are found by differentiating equation ilOS ) . 


V(t) = Vo + Vot + 7 + . • . 


B^U 2 

(110a) 

> 

3 

1! 

O 

> 

Vo = Vo = Vo = 0 

(110b) 

iv _ 2 b2o) y 

Vo - ^0^0 

(110c) 


Next, an interaction potential U*Ct) will be found which simultaneously 
produces the correct vibrational excitation and satisfies the conservation of 
total energy 


E = T + U* + V 


Again the assumption is made that the potential is of simple exponential form 

U*(t) = U* exp |- -]| (111) 

but r(t) in this equation is distorted from its usual adiabatic collision 
value. Then 

u* = _ £ u* (112) 


which may be combined with 


T = 


mf^ 

2 


to give 



(113) 


This expression for the kinetic energy in the conservation equation yields 
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E 


(114) 


V(t) = 



This equation is evaluated at the turning point and combined with equa- 
tion (110a) to give 


Vo = E - 



(115) 


But equation (112) shows U* = 0 since 1 -^ = 0, and so equation (115) may be 
solved for U*. 



(116) 


This same procedure may be used to find values for higher derivatives of U* 
at the turning point when equation (114) has been differentiated repeatedly. 
For example, after it is differentiated twice. 


V = mL^ 
o 


5U*US 

/ii$f 


/UoV 

L^o 

" \ u* i 
\ 0/ 

u* 

o 



- u; 


(117) 


This equation is solved for U* by using 


Vo = 


0 from equation (110b) . 


Similarly, 


Uo = 0 

,,,iv 4US^ 2 


(118) 


(119) 

( 120 ) 


and higher derivatives can be found, if necessary. 

Return, now, to the exponential dependence of U*(t) on r(t) in 
equation (111). For an adiabatic collision trajectory, 

j I = sech^ at 

with a = u/2L. This makes plausible the functional form for nonadiabatic 
collision, that is, one where the internal modes of the colliding partners 
take part in the energy transfer 
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|- " 4>(t)sech2 


at 


when r(t) is distorted by transiently stored vibrational energy. Therefore, 
it is assumed that 


U^(t) = Uo(j) (t) sech^ at 


( 121 ) 


The function 4)(t) is near unity since the collision energy is large when 
compared with the stored vibrational energy, and it is not unduly restricted 
by the assumption, physically realistic, that ())(t) ^ constant as 1 1 1 ^ 

This last assumption is used with the derivative of equation (121) to give 
a = u/2L as before. 

The transition probability, P*, is influenced most strongly by the 
behavior of U*(t) near the turning point, and the expansion 


U*(t) = Uq sech^ at Uq + (p^t + j + 


1 


( 122 ) 


will be sought. The first coefficient is found to be 


= 1 


(123) 


by evaluating equation (121) at the turning point. The first derivative of 
equation (121), also evaluated at the turning point, reduces to 




and since U = 0, it follows that 
o 


= 0 


(124) 


Further differentiation of equation (121) and use of equations (118) , (119) , 
and (120) gives the expressions 


n* II* 

= uf ^ ^ . 2a^ 


= 0 


iv ^ 2 

° ■ m2L‘* ■ 3 hmL2 


+ 12a2*’ - 16a^ 


(125a) 

(125b) 

(125c) 
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The expansion for can be written 


<j)(t) = 1 + c(at)^ + dCat)*^ + . . . 


where 


and 


c 


<i>0 

2I2 


d = 


24a‘* 


and will be terminated after the term containing t‘^. This is justified by 
the fact that only the behavior of ((iCtj inside the interval | at | <_ 2 is 
important. Outside this interval, sech^at changes very slowly, and since 
<|)(t)sech2 at must remain close to sech^ at, its Fourier transform 


F* 



(j) (t) sech^ (at) e^^^^ 


dt 


(126) 


receives nearly its full magnitude from -2 ^ at <. 2. Therefore, it is taken 
that 


U*(t) = Uq sech^(at)[l + c(at)^ + d(at)‘*] (127) 

is the interaction potential which produces the correct amount of stored 
vibrational energy and, at the same time, satisfies total energy conservation. 
The evaluation of its Fourier transform, F*, to obtain the transition 
probability proceeds by contour integration. 

In order that F* may be evaluated, the integrals 

f = r t^ sech^(at)e^^^ dt = -^ T sech^ dx (128) 

•'-00 •'-00 ^ ' 

must be found. But these integrals can be treated in the same manner, using 
the same contour and limiting process, as is used in appendix B for a slightly 
different integrand. The residue at the isolated singularity inside the 
contour is 

(If 

and for ca/a >> 1, 
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so that 


n+ 1 . n 

^ o -n ^ ^ ^ - TTa )/ 2 a 

fn = 2TriRjj » e 


^n-i n+2 
2 a 


- 2. (I) 


(130) 


and 


= 4tt^ 


(I) (i - T - #) 


If P* is used to denote the corrected transition probability, after the 
temporarily stored vibrational energy is accounted for, then 


P* = 3^ , fp-') |f*| 

w Vhco/ ' ' 


(131) 


and 


P lF|^ 

However, when equation (130) is compared to equation (46) whose right hand 
side is F/to, this ratio is seen to be 



^ TT^C 

Tr‘*d\^ 

^ ■ V 

^ ■ 4 ■ 

16 j 


and therefore 




where P is the transition probability of equation (47), and where 


c = 


d = 


hto 
2 


-1 + 


©■]-i (¥)©(;) 


- + C 


The effect of the correction factor is to make, roughly. 


(132) 
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(133) 


p* 



Consequently, the collision cross section, S* (eq. (4)), and the rate 
coefficient, a* (eq. (5)), are diminished by this same factor. But, the 
decreased rate coefficient implies (cf. eq. (6)) a relaxation time, x*, which 
is about 5/2 times x, the relaxation time predicted by the theory which does 
not incorporate the effect of the temporarily stored vibrational energy. 


COMPARISON WITH EXPERIMENT 


There is no lack of data on rates of vibrational relaxation as a function 
of temperature (ref. 7). The experimentally determined quantity is the relax- 
ation time X, and according to equation (6) the product of x and the 
pressure p is 

kT 

‘ a.„[l - 

In this equation ajo is the rate coefficient for transitions from the first 
excited state to the ground state, given by equation (95). This connection 
between the analytically derived a and the experimentally obtained t pro- 
vides a means for testing the analysis. At the same time, it yields a method 
for determining effective values for the molecular potential parameters, A 
and L. 

The procedure begins with selectin_g a molecular species and collision 
partner and assuming trial values for A and L to describe the intermolecular 
potential. Thereupon, the rate coefficient a^oCT, A, L) from equation (95) 
determines, through equation (100), a relaxation time x which can be com- 
puted for chosen values of temperature and unit pressure. These computed 
values will be (nearly) collinear on a Landau-Teller plot, and the slope_^and 
ordinate of this line can be adjusted by changing the trial values for A and 
L. A combination of values for A and L can be determined by eye which pro- 
duces the best agreement between the theoretical line and experimental data. 
This combination is taken to be an acceptable approximation to the effective 
values of these parameters. An example of the result of this procedure is 
shown in figure 6 for 02 -Ar collisions and in figure 7 for N 2 -N 2 collisions. 

It will be noticed in figure 7 that the slope of the px line is not 
quite correct for perfect agreement between theory and experiment. This 
observation leads to the suspicion that further adjustment of A and L will 
yield even better agreement than is shown, but this is not the case. Further 
adjustment to improve the slope moves the line vertically and degrades the 
overall agreement. Also shown in figure 7 is the relaxation time predicted 
for the delta function potential calculated by Meador (ref. 30) . This 
predicted relaxation time is the order of 100 times larger than observed, and 
some remarks will be made later about this difference. 
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Figure 6.- Landau- Teller plot of vibrational Figure 7.- Landau-Teller plot of relaxation time 

relaxation for 02~Ar collisions. for N2-N2 collisions. 

The above procedure has been applied to a variety of diatomic species 
undergoing collisions with an atomic or diatomic collision partner. The 
results are given in table l_where two sets of intermolecular potential param- 
eters are shown. The pair A,L were deduced using the semiclassical method 
without correction for conserva^tion of total energy, and when this correction 
was made, their counterparts A*,L* were obtained. Comparison of these pairs 
shows that the sizable difference found earlier 



produces only slightly deviating forms for the exponential intermolecular 
potential. Close numerical agreement of the two sets of constants makes it 
difficult to ascertain which pair furnishes better agreement between experi- 
ment and theory. Hence, the pair A,L produced by the more simple, 
uncorrected semiclassical theory will be used in the remainder of this paper. 
Care must be taken in using these values for the constants at extrapolated 
temperatures. Neither the experimental data nor assuii5)tions made in the 
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analysis can be meaningfully extended to temperatures sufficiently high for 
dissociation to become a competing process with vibrational excitation. 

The accuracy of the values given for A and L cannot be assessed 
quantitatively_^ A possible approach to constructing table 1 is to find com- 
binations of A and L that make the px line fit the experimental data best 
in the least squares sense. However, the numerical problem posed by this 
approach has proven to be without a distinct solution. There is a spectrum of 
combinations of A and L that appear to be equally valid from a numerical 
least squares standpoint. Other numerical attempts have also proved 
disappointing, and it was necessary that the px line be fit to ttae data 
subjectively. However, experience in choosing trial values for A and L 
gives assurance that th^ '^best" value for L can be determined within 
10-25 percent and for A within a factor of 2 or 3. In general, the more 
stringent confidence limits apply for heavier collision partners. 

Attempts to evaluate A and L for the halogen gases were without reward 
for Br 2 ~ Br 2 and I 2 - I 2 collisions, and there is near-indeterminacy for 
C^2 " C^2* "The experimental data C^^f- 31) for these gases are obtained by a 
measurement of sound absorption at temperatures very low relative to the other 
data. A Landau-Teller plot of the halogen data is, an£malously, a nearly hori- 
zontal line. Except for chlorine, no combination of A and L predicts a 
slope small enough to give meaningful agreement, and even for CI 2 the values 
of L required to fit the data seem to be abnormally small (the order of 
0.05 A) . 

A direct numerical comparison of experimental and calculated px-values 

is m£de in table 2. The computed values were obtained using the combination 

of A and L found in table 1. It is to be noted that the comparison is made 

at temperatures which lie at or near the ends of the temperature range where 

A and L are fit to the experimental data. These are the temperatures where 

px and px differ most (see figs. 6 and 7) in that range. Conse- 

e xp XX c a X c 

quently, table 2 gives bounds for the discrepancy between theory and experi- 
ment. Also, since px varies linearly with on a Landau-Teller plot, 

table 2 can be used to construct a counterpart to figures 6 and 7 for any of 
the species listed. 

Strictly speaking, the CO data should be evaluated with a more general 
theory for heteronuclear molecules, since the effective potentials for C and 
0 atom interactions can be different. Corrections for the asymmetrical cen- 
ter of mass should be small in this case, however, and the homonuclear mole- 
cule theory may be expected to give an approximate average interaction 
potential, corresponding to the random orientation of the molecular axis that 
occurs in the experiment. Dipole interactions affect the long range forces, 
but probably do not contribute largely to the steep part of the potential 
responsible for vibrational transitions. 

Fitting vibrational relaxation data is widely conceded to be a rather 
inaccurate way to determine intermolecular potentials. Nevertheless, such 
potentials are so urgently needed for analysis of a wide variety of atomic and 
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TABLE 1.- VALUES OF THE POTENTIAL PARAMETERS (A,L) AND (A*,L*) 


Species 


O2 

O2 

O2 

02 

N2 

CO 

CO 

CO 


02 

N 2 


TABLE 2.- 


Species 


O 2 


O 2 

O 2 

O 2 

N 2 

CO 

CO 

CO 


Partner 

A(eV) 

L(A) 

A* (eV) 

L*(A) 

Approximate range 
of T, °K 

O 2 

800 

0.19 

900 

0.18 

500-5,000 

Ar 

800 

.24 

900 

.23 

1,000-5,000 

He 

550 

.26 

600 

.26 

300-2,000 

»2 

300 

.27 

300 

.26 

300-2,000 

N 2 

700 

.23 

800 

.21 

2,000-8,000 

CO 

850 

.19 

900 

.18 

1,200-6,000 

H 2 

300 

.28 

350 

.27 

300-2,000 

He 

400 

.27 

500 

.26 

300-2,000 


Calculated delta function model 


O 2 

147 

0.47 



1,000-10,000 

N 2 

350 

.40 



1,000-10,000 


COMPARISON OF EXPERIMENTAL AND CALCULATED RELAXATION TIMES 


Collision 

T, °K 

-T»- 1/3 

pr , atm- 

-sec 

partner 

1 

Experimental 

Calculated 

O 2 

5000 

0.058 

4.0x10-7 

5.8x10-7 


2000 

.079 

6.0x10-6 

7.2x10"6 


1000 

.1 

1.0xl0-‘t 

1.2x10-*+ 


500 

.126 

2.5x10-3 

2.2x10-3 

Ar 

2000 

.079 

3.0x10-6 

3.1x10-6 


1000 

.1 

8.9x10-^+ 

1.2x10 3 

He 

1500 

.087 

1.0x10-6 

0.5x10 6 


300 

.149 

4.0x10-6 

6.0x10-6 

H 2 

1000 

. 1 

3.8x10-7 

3.2x10-7 


300 

.149 

2.4x10'6 

3.7x10-6 

N 2 

5000 

.058 

6.6x10-6 

9.3x10"6 


2000 

.079 

7.4x10-*+ 

5.3x10"*+ 

CO 

5000 

.058 

1.9x10-6 

2.5x10"6 


2000 

.079 

6.8x10-6 

6.6x10~6 

He 

1000 

.1 

2.9x10-6 

2.7x10"6 


300 

.149 

2.6x10-3 

3.0x10-3 

H 2 

1000 

.1 

2.8x10"6 

2.9x10"6 


300 

.149 

6.8x10-6 

6.8x10"6 

. 
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molecular problems, and so little information about these potentials exists, 
that there is always a strong temptation to evaluate the potentials as done in 
table 1, just to see what relations appear. 

First of all, the results are generally consistent with similar attempts 
made with one -dimensional vibrational excitation theories Any one- 
dimensional theory can only be used to determine L, not A, since it can fit 
only the slope and not the absolute magnitude of the data on a Landau-Teller 
plot. Schwartz, Slawsky, and Herzfeld (ref. 11) obtained reeisonable agreement 
with data (within a factor of 20 or so) when they chose L to fit the 

repulsive part of a Lennard-Jones potential, U = e[(rQ/r)^^ - Th^y 

found that L r^/lT.S, and for values of Tq determined from relatively low- 
temperature viscosity data (ref. 32) this gives L - 0.2 A for O 2 , N 2 , and CO. 
Benson and Berend (ref. 33) chose to fit the repulsive part of a Morse poten- 
tial, U = Uo[(l - e'^^)^ - 1], and they presented values of a that best fit 

measured vibrational transition rates. In this case, L - l/2a, or L - 0.25 A 
for O 2 , N 2 , and CO. Rapp and Sharp (refs. 21,22) and Treanor (refs. 25,26) 
both used L = 0.2 A for their calculations of N 2 -N 2 vibrational excitation. 
These values are about the same as given in table 1 and are within the 
uncertainty inherent in fitting the data. 

However, it seems to have generally been overlooked that all of these 
results from vibrational relaxation theory are different from estimates of 
N 2 -N 2 and O 2 -O 2 potentials that have been deduced by other methods. Mason and 
Vanderslice (refs. 34, 35) developed a delta function model for calculating 
short range intermolecular forces, and Meador (ref. 30) refined this method 
and applied it to N 2 -N 2 collisions, showing that it agrees well with the 
potential for such collisions given by Amdur, Mason, and Jordan (ref. 36) , who 
deduced this potential somewhat indirectly from Ar-N 2 and He-N 2 scattering 
measurements. Meador (ref. 30) then applied this same delta function model to 
O 2 -O 2 collisions; in this case experimental verification is lacking. Meador’s 
values of A and L are listed in table 1, and it can be seen tha^ L is con- 
siderably larger than deduced from vibrational relaxation while A is smaller. 
At the collision energies of interest (0.1 - 1.0 eV) , the delta function poten- 
tial is less steep and of longer range. The shallower potential results in a 
more adiabatic collision, smaller excitation probability, and larger relaxa- 
tion time (as shown in fig. 7). The expansions involved become doubtful for 
these large values of L, so the exact numerical results for Meador’s poten- 
tial should not be taken too seriously, but the qualitative effects should be 
valid, at least. 

The disagreement between the potentials is not surprising, when one 
considers that a multiplicity of interaction potentials exist for molecular 
collisions, just as for atomic collisions, depending upon how the orbital and 
spin momentum vectors add up. Scattering is primarily produced by the long- 
range outer potentials that correspond to the largest total spin and degener- 
acy, such as calculated by Meador (ref. 30). On the other hand, vibrational 
transitions are undoubtedly due to interactions that proceed along shorter 
range, steeper potentials that occur when some of the spins are paired, and 
vibrational relaxation rates are presumably telling us something about the 
shape of these inner potentials. However, a number of factors prevent us from 
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deducing precise quantitative values for these inner potentials from 
vibrational relaxation data. Perhaps the most serious is that the semiclassi- 
cal method uses the spherical average of the true potential and does not 
account for energy adiabatically transferred to and from the vibrational mode 
during the course of the collision. This effect is dramatically shown in the 
numerical calculations of Sharp and Rapp (ref. 22) . The perturbation is not 
really the simple adiabatic scattering impulse assumed in equation (35) , even 
though this shape has been traditionally used in all semiclassical analysis of 
the problem, but is distorted by the temporary transfer of energy into high- 
order vibrational modes, and is also slightly skewed by that small portion of 
the energy transferred irreversibly. Even where changes in impulse shape are 
relatively small, the effects on the Fourier transform can be sizable. Thus, 
the potential deduced from vibrational relaxation measurements should be 
thought of as an "effective” potential which incorporates the effect of energy 
stored in the vibrational mode, and it should not be interpreted as a pure 
interaction potential. 

Several other approximations involved in the theory may be sources of 
error also, particularly at higher temperatures. One of these is due to the 
fact that small perturbation probabilities become the order of 0.1 at the 
velocities for the highest temperature data, and at this point become 

sizably greater than the actual transition probabilities (refs. 24, 26). How- 
ever, this error is compensated because multiple step transitions occur at 
these temperatures, and the total energy transferred is the same in either 
case. A second error which decreases the relaxation rate at high temperature 
occurs because the upper limit of the integral for the cross section is then 
finite, leading to the result of equation (60), rather than the approximation 
of equation (61) which is valid for very small transition probabilities. A 
third error occurs at high collision energies because the expression for the 
perturbation potential (eq. (19)) should include higher order terms in the 
expansion of the distances ri and r 2 . Finally, the real interaction poten- 
tial_^ are not simple exponentials. Thus, different effective values of L 
and A could be used for different ranges of temperature and adjusted to fit 
data exactly. 

In spite of the uncertain interpretation of the potentials used, the 
simple theory can be fit to observed vibrational relaxation data reasonably 
well. The best fit is excellent when the colliding partner is an inert atom 
such as Ar or He, which is the type of collision that the theory is designed 
to describe. The correlations are only a little worse for collisions such as 
O 2 -O 2 and N 2 -N 2 . The principal point to be emphasized is that the "effective 
potentials" deduced from these correlations are not necessarily realistic for 
other purposes. 


DISCUSSION OF RESULTS 


Perhaps the most satisfying result from this study of vibrational 
excitation is that the problem has been treated analytically in three dimen- 
sions, and the cross sections and the rate coefficients are expressed as 
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relatively simple functions of temperature and the interaction potential 
parameters only. To first order, the results show the same fionctional 
behavior as one-dimensional models. This behavior confirms the qualitative 
results in reference 16 and further justifies the widespread use of the 
one-dimensional model for analyzing vibrational relaxation data. 

The relative simplicity of the semiclassical model compared with the 
method of partial waves (which gives slowly converging results in three- 
dimensional quantum treatments) has made it possible to retain fully analytic 
expressions all the way to the final rate coefficients. The results clearly 
show how vibrational transitions (predominantly Av = ±1) are accompanied by 
simultaneous rotational transitions (predominantly AZ = 0, ±2), an aspect of 
the process that has often been ignored. The effect of these coupled rota- 
tional transitions increases the total vibrational transition rate from 50 to 
100 percent at typical conditions. This effect went unnoticed in one- 
dimensional theory because it could simply be absorbed in the empirical steric 
factor involved in that theory. Herzfeld and Litovitz (ref. 37) and others 
have pointed out that such coupled rotations must occur in three-dimensional 
collisions, however. Also, the three-dimensiona^ model permits one to evalu- 
ate effective values for both the scale factor A and the range factor L of 
the intermolecular potential by fitting vibrational relaxation data, whereas 
one-dimensional models permit only the factor L to be determined. Both _ 
models give about the same values for L, the order of 0.2 A. The value of A 
given by the three-dimensional model is reasonable, the order of several hun- 
dred eV for typical interactions. These potentials are much steeper and of 
shorter range than known scattering potentials, and we suggest that this 
occurs primarily because vibrational transitions are produced most frequently 
when the interaction proceeds along one of the shorter range, steeper poten- 
tials that obtain when some of the electron spins involved are paired. How- 
ever, the ^effective" potential for vibrational relaxation should probably not 
be interpreted as a pure interaction potential in any case, since it also 
includes the effects of potential stored in high vibrational modes during 
collision. 

Although the interpretation of the potentials derived may be questionable, 
the semiclassical theory of vibrational excitation can be fit reasonably well 
to most relaxation- rate data if the appropriate effective potential is chosen. 
The theory appears to contain enough realism to be suitable for investigating 
some other facets of the vibrational excitation problem; for example, the 
effects of anharmonicity , vibration-vibration exchange, heteronuclear 
molecules, and perhaps triatomic molecules. 

The square of the transition matrix elements decreases in size fairly 
rapidly as |aZ| increases, but the transition probability increases rapidly 
if Av and AZ have opposed algebraic sign because the energy difference 
between states involved in the transition is then decreased. The combined 
effect of these two factors is that transition probabilities for |aZ| = 0,2 
are about equal, while those for |aZ| = 4 are about an order of magnitude 
smaller. Higher order perturbation terms can increase the later transitions 
somewhat and contribute some higher order, even number changes AZ as well. 


40 



but it appears that AZ = 0, ±2 will remain the dominant rotational 
transitions which are coupled to the vibrational transitions Av = ±1. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, April 8, 1970 
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APPENDIX A 

EVALUATION OF THE CONSTANTS 3^(Av, hV) 


In this appendix, the overlap integral 


6(v,Z,m;V ,1' ,m') 



(V 


f 62 /P-Pe\/ . 2o 2 

' ,Z' ,m')Lsinh 6 \ 2L ® 


<!> 


sin^ 0 cos^ (f) 


')] 


^v,l,m)P" ® 


(Al) 


will be evaluated. The square, (v, Z.,m; v' , Z’ ,m' ) e B^ (v, Z,m;v+Av, Z.+AZ,m+Am) 

will then be averaged over all possible values of m and, with the assumption 
that 1 >> 1, it will be found that this average is independent of Z. Thus, 
for vibrational transitions from a fixed vibrational level, the averaged 
squares are functions of Av, AZ-, and Am only; B^(Av, AZ, Am). Since these 
squares are proportional to the transition probability, they may be combined 
in the manner probabilities are combined. 


B2(Av, AZ) = b'^(Av, At, Am) 

Am 


EVALUATION OF THE OVERLAP INTEGRAL 


The unperturbed wave functions ^ j^^CPjBj 4>) and ^ (p 

that appear in equation (Al) are, respectively, the wave functions of the 
rotating and vibrating diatomic molecule before and after the collision. For 
a rigid- rotator, harmonic oscillator these wave functions can be written: 

( n -(a/2)Cp-Pg)^ 

(v, Z ,m) ^ Hy[/a(p-Pg)] 


In this equation, Hy[/cT(p-P^) ] is the vth Hermite polynomial and the 
Y^^(6,(|)) are spherical harmonics expressed in terms of (cos 0), the asso- 

ciated Legendre functions of the first kind. 







imi 

Z 



(cos 


S) 


The normalizing constants are expressed 


42 



iml 


f2Z+l {I - imon 

L 2 (^ + imi) ! J 


1/2 


(A2) 




1/2 


(A3) 


where 


a = 


yco 


It is seen that equation (Al) can be written 

B(v,Z,m;v' ,Z' ,m') = Ii I | Y^,^, (6, 4 >) (sin^ e 

•'6=0 •'4>=o '• 

r 2 \ 

+ -g- sin^ 0 cos^ Y2^^(0,(j))d0 d(f) 


cos" 


(A4) 


with 


Ii = 


2La 


^ CyCy I I (x) Hy (x) xe dx 


(X = v/^(p - p^)) 


The orthonormal properties and recurrence relation for the Hermite polynomials 
which will be used to evaluate Ix can be found in mathematical texts (e.g., 
ref. 38, p. 308). These properties are 


Hv* (x)Hv(x)e'' dx = 0 (v / v’) 

= 2^v! /IT = /a/Cy^ (v = v') 


(A5) 


and 


xHy(x) = ^H^^.^(x) + vHy^x(^) 

.2 


(A6) 


When this last equation is multiplied by e (x) and then integrated over 

( 0 , <») , there is obtained 

X oo 2 X 2 /•“ 2 

Hv'(x)Hy(x)e ^ X dx = - J Hy+j (x)Hv.(x)e”^ dx + vj ^ (x)e^ dx 


(A7) 
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When equation (A5) is used with equation (A7) , it is seen that the right-hand 
side of CA7) is zero unless v' = v ± 1. Therefore, 


C ,C 
v' V 

■ ^ 2La 


/v + 1 

/ 2a 

when 

v'=v+l 

(A8a) 

/ 2a 

when 

r-l 

1 

> 

If 

> 

(A8b) 


when 

v' ^v±l 

(A8c) 


The integral of equation (Al) may now be written 




(A9) 


where (v + 1,0) is equal to v + 1 or v accordingly as v’ = v + 1 or v - 1, 
and where 


^ 2 7T 

= ^ I 


^ ^ cos^ cj) dc() 


plm* I 1ml 




(cos 0)sin^6P^^ (cos 0)d0 


(AlO) 


and 


6 ^ fl i(m-m')(() 4 , ,,1 I m' I „ I ml qa • Snnin'l/- 

= -g- 2^ f e ^ ^ cos^ <t> dcj) C^, I P^, (cos 9 )sin^e (cos 0 )d 6 

*^o J L •'o 


(All) 


The method for evaluating each of these integrals will be demonstrated 
separately. In order to facilitate notation in what follows, |m| will be 
written as m. 


(i) I 2 


^ 2tt 

1 I i 

2tt I ® 

- •'o 


(m-m')d) ^^^2 


cos^ dcf) 


..m 


'em r pm' 

'^7 I ^7' 

^ r\ 


o m 

(cos 6)sin^ 6 P^(cos 0)d0 


The first integral in equation (AlO) can be written 
2tt . 


1 i 

2tt J ® 

•'o 


^ cos^ (j) d(j) = - — 


2tt 


X 7 T ^TT 

cos (m-m’)(j) cos^ cj) d(f) + i I sin(m-mOc|) cos^ cf) d^ 

7T •'-TT 
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The second integral on the right-hand side has symmetric limits and an 
integrand which is an odd function. Its value is therefore zero. To evaluate 
the first integral write cos^ ((> = (1 + cos 2<j))/2. Then 


1 i(m-m*)(|) 

27T J ^ 

•'o 


COS^<fi 



cos(m-m’)(j) d(|) + ^ 



cos 


2(|> cos(m-m*)(l) d<i> 


The first integral on the right-hand side is equal to zero unless m’ = m, in 
which case its value is 27 t. The second integral is zero unless m’ = m ±2, 
when its value is tt. Hence, the chosen form of the interaction potential 
(eq. (19)) permits only the changes m’ = m, m + 2, or m - 2 by a collision. 
Accordingly, 


^ 1 
2 

I 

4 


if Am 
if Am 


0 

±2 


pm* m 




0)P7(COS 


0)sin^ 0 d0 


If X = cos 0, then 


,C^ f P^, (cos 0)P^(cos 0)sin^ 0 d0 = I 2 ^ 6 = C^, C™ f P^, (x) (1 - x P™(x) dx 

•'n - 1 


'O - 1 

The associated Legendre functions are orthogonal in (-1, 1) : 


p“, (x)p^(x)dx = 0 a' ^ 1) 

= (C™)"2 (^. = 

They also satisfy the recurrence relations 

/x^ - 1 P™^^ (x) = il - m)xP™(x) - (Z + m)P™ ^ (x) 


(A12a) 

(A12b) 


(A13) 


xP^(x) 


1 

2Z+1 


(Z- m + l)P™^^(x) + (Z + m)P™_^(x)j 


(A14) 


P®^l(x) = P“_^U) + (2Z+1)4^-1 P™"^(x) 


(A15) 


These rel ations h old for real and complex x. For x in [-1, 1], the 
radical /x^ - 1 which appears in (A13) and (A15) is complex. However, these 
relations will be used to develop formulas involving (1 - x^) which is real in 
the interval. 
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It has been found that Am can take on three possible values in I 2 , 
(Am = 0, ±2), and for each of them, I 2 q will have a separate value: 

(a) I q(Z,Z';0) = C®, cj f P™, (x)(l - x2)p^(x)dx 

> •'-1 

(b) I^^q(Z,^'; + 2) = P^;^(xKl - x2)p“(x)dx 

(c) I^ q(Z,Z';- 2) = P^;^(x)(l - x2)P™(x)dx 

’ *'-1 

Each case must be considered separately, and notation will further be 
abbreviated by writing P^(x) = P?. 


case (a): I,_e ' ^".C^ 


7_‘ - /I 


dx 


When the orthogonality property (A12) is applied to the first integral on 
the right-hand side of the above equation, it is seen that this integral is 
zero unless - I when its value is 1/(C^)^. 

The second integral on the right can be evaluated by multiplying the 
recurrence relation (A14) by (xP^^), and then applying this same equation to 

the terms 1 which appear in the result. Finally, it is 

obtained that 




,m 




Z + m 
2^+1 \2Ul 


1 


[ (1 - m)P^P™, + (Z. + m 


1)P^,P^ 


- 2^1 


Z--m+l 

2Z+1 


1 

2Z+3 




(Z + m+ 1)P",P^ ) 


(A16) 


When both sides of this equation are integrated over (-1, 1), the orthogonal- 
ity relation (A12) shows the right side of the resulting equation to be zero, 
except when V =ZorZ+2orZ-2; that is, for aZ = 0 or ±2. The 
results are: 


q(AZ= 0, Am=0) 


, (Z+m) (Z-m) (Z-m+1) (Z+m+1) 

(2Z+1) (2Z-1)" (2Z+l)(2Z+33 


(A17a) 
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Case (b) : 


■-[ 

■ {“: 


(Z-m+1) (Z-m+2) (Z-+m+2) (Z+m+1) 


(2^+1) (2J+3)2(2Z+5) 




m) (Z^+m-1) (Z-m) (Z-m- 


(2^+1) (2Z-1)2(2^-3) 




1/2 






(1 - 3C^)P^ dx 


(A17b) 


(A17c) 


In this case, the recurrence relations (A13) and (A14) are used to write 
P^, in terms of P 2.»+2 ^7'-2‘ preliminary result of combining 

these two equations is 

- /x2 - 1 [(Z + m+ l)(7 + m)P^_^ - (Z- m)C^- m + DP^^J (A18) 

and from this equation comes 

- /x^ - 1 P^-"^ = [(Z + m+ 2)(Z + m+ Dp’^^J - (7- m - 1) (Z - m)P^;^ j] 

This last equation is then mu ltipli ed by /x^ - 1 a nd equation (A18) is 
applied to terms involving /x^ - 1 p’?'*’^ and /x^ - 1 Finally, it is 

derived that 

ri _ _ (Z'+m+2) (Z'+m+l) f (Z'+m-l) (Z'+m) ,^m (Z' -m-1) CZ^' -m) ,^m H 

Z’ " 2Z'+1 L 2Z'^1 ^Z'-2‘ 2Z’-1 PZ'J 

(Z ' -m- 1) (Z* -m) r (Z ' +m+l) (Z ' +m+2) m (Z'-m+l) (Z'-m+2) ^^m ”| 

2Z'+1 [ ‘ 2Z'+3 ^Z'" 2Z'T3 ^Z'+2j 

(A19) 

This equation is multiplied by P^ and the resulting equation is integrated 

over (-1, 1). It is seen once more that the orthogonality relation (A12) 
precludes rotational transitions in which Z' is not equal to Z or to 
Z ±2. The results are: 


= -2 


(Z+m+2) (Z+m+1) (Z-m) (Z-m- 


(2Z+3)2(2Z-1)2 


if 


1 /2 


(A20a) 
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(A20b) 


Case (c) : g 


(Z+m+4) (l+m+3) (Z.+m+2) (Z-+m+l) 
(2^+5) C2Z+ 3)2(21+1) 


1/2 


(Z-m-3) (Z-m-2) (Z-m-1) 
(2Z+3) (2Z-1)2(2Z 


)C^-m)' 

+ 1 ) 


1/2 


{I' ,l\-T) = C^.^C™ 


•'-1 


- x2)p^ 


dx 


(A20c) 


Here, the recurrence formula (A15) is used repeatedly to obtain 


-F^U ^ - C2l’ + 1) Z I* 


97t-l ^7^7T 


1 pmpm 


+ 2~2Z'+3 ZZ' 2Z'-1 "Z"Z' 2Z’-1 Z"Z 


■4 


and this equation is then integrated over (-1,1) to get 


^2 ^ 


(Z+m) (Z+m-1) (Z-m+2) (Z-m+1) 


"1 1 /2 


(2Z+3) (2Z-1) 


(A21a) 




(Z-m+4) (Z-m+3) (Z-m+2) (Z-m+1) 


-. 1/2 


(2Z+1) (2Z+3)2(2Z+5) 


(A21b) 


I ,^ g (- 2 ,- 2 ) = - 


(Z+m) (Z+m-1) (Z+m-2) (Z+m-3) 


ni/2 


(2Z+1) (2Z-1)2(2Z-3) 


(A21c) 


This completes the evaluation of I 2 (Z, Z' ;m,m' ) = l 2 (AZ,Am), and it is 
seen that this integral can take on a total of nine different expressions, one 
for each combination of the three allowed changes, AZ, in the rotational 
quantum number and the three allowed changes. Am, in the magnetic quantum num- 
ber. If it is assumed that Z >> 1, each of these expressions can be simpli- 
fied. For example, a reasonable approximation for ^(0,0), using Z >> 1, 
is 

(Z+m)(Z-m) (Z-m+1) (Z+m-2) 2(Z2-m2) _ Z2+m^ 

■ (2Z+1)(2Z-1) ' (2Z+l)(2Z+3) ^ ~ /[il 2l^ 

Therefore , 

l2(0.0) = i I2_J,C0,0) (A22a) 
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Similarly, 


12 ( 0 , 2 ) 


l2(0,-2) 


12 ( 2 , 0 ) 

12 ( 2 , 2 ) 


l2(2,-2) 


l2(-2,0) 


l2(-2,2) 


l2(-2,-2) 


1 

4 

1 

4 

1 

2 

1 

4 

1 _ 

4 

1 

2 

1 

4 

1 

4 


2,0 


2,0 


■2.9 


2,0 


2,0 


2 . 0 


■ 2,0 


2,0 


( 0 , 2 ) « - 


( 0 ,- 2 ) 


Z^-m 2 

8^2 

i2-m2 

8^2 


( 2 , 0 ) 

( 2 , 2 ) « 


Z2-m2 

8^2 

{Z+m} 2 

16^2 


( 2 ,- 2 ) 


16^2 


(- 2 , 0 ) 


Z 2 -m^ 

8^2 


(- 2 , 2 ) 


1612 


(- 2 ,- 2 ) 


(l+m) 2 
16^2 


(A22b) 

(A22c) 

(A22d) 

(A22e) 

(A22f) 

(A22g) 

(A22h) 

(A22i) 


r • • 'v T 
(XI) l3 = -^ 


2tt 


X 


27T 


cos 


i (m-m* ) (j) 


d(j) 


^m* m 


f 

n 




(cos 


0) (sin^ 0)P^(cos 0)d0| 


The first integral on the right in the equation for 1 3 can be written, 
with the use of cos^ tf) = (l/ 8 )cos 4 (}) + cos^ cfi - 1 / 8 , as 


J_ 

2 tt 



4 

COS^ 


4, 



COS^ 


(p cos(m-m*)(j) dcf) 



cos (m-m’)(|) d(p 


+ 



cos 


4(p 


cos (m-m* J(p d(p 


It was found earlier that 


2 tt 



cos^ 


cos (m-m’ ) (j) d(p 



if m* = m 
if m* = m ±2 
otherwise 
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Hence, after evaluating the second integral 


if. 


2 TT 


if m» = m 

cos** d(f> = if m' = m ±2} + 

0 otherwise 


j_ r 

16tt J 

•'-ir 


cos 4<fi cos Cm-m* )(j) d<j> 


The last integral vanishes unless (m - m*) =4, in which case its value is 
equal to ir. Consequently, 


«2ir 


*'o 


i (m-m ' ) 4> 


d<t> = T 


3 
8 
1 

4 
1 

16 


if 

m* 

= m 


if 

m* 

= m 

±2 

if 

m’ 

= m 

±4 

m' . 





(A23) 


When X = cos 0, the integral 


I 3 e=c^;c^ r P^lccos 0 )(sin 5 e)P^(cos 0)de=c™;c“ f' P^| (x) (1 - x^)2p^(x)dx 
’ *^o •'-1 


Now 


(1 - x2)2 = (1 - X^) - X^(l - X^) 


so that 1^,0 becomes 
3,0 




dx 


= I 


:,6-'=Z-‘=z/’ P”![^^(1-*^)]P” ax 


Since I,. _ has already been evaluated, it remains only to evaluate the 
integral 


'j.e = - x2)lp” dx 


CA24) 
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In order to carry out this integration, the recursion relations obtained from 
equations (A15) , (A16) , and (A19) will be used. For example, if we use 
Z » 1, equation (A16) may be written 


[2Cl2-m^)P^,p”* (l*m)2p”,P^_^. (I-n02p^,p“^J (A25) 

This equation is multiplied by to give 

^ [2(Z2_ n,2)P™,x2p^+ + ^^^6) 


In similar ways there are obtained 


P^yx2(l-x2)p”^ jij. [(Z*m)“p5',_^x2p^-2(Z2-n,2)2p^,x2p“. (Z - m) -p” ^^2p^ ] 


P"’x2fl - x2-)p"'-2 r2p"'x2p"' -p"’,(2p"' - p"’v2p™ 

I ~ \r^l ^l' Pz^P^»+2 P^^ PZ'-2_ 


(A27) 


(A28) 


P^(l - x^)2p"'; 


^ 4^2 


^ [2P“(1 - x2)P®;2 - p^(l _ x2)P^;2^ - p®(l - x2)P^;2^j (A29) 


It has been found that Am can take on five possible values in 1 3, 
(Am = 0, ±2, ±4) and for each of these, I3 q will have a separate value: 


Case (a): I ■ (AZ,0) = C^, C™ j‘ P™,x2p” dx-C™,C^ pj, 


.x'tp'? dx 


When equations (A25) and (A26) are integrated over (-1,1) and the usual 
orthogonality conditions are applied. 




^2 - m2 (j2 - m2) 2 f2’2 - m2)' 


2J2 






Therefore 




Z-2-m2 Z.2-m2 ^ ^ (Z2-m2)2 

2Z2 " 2^2 8 ^4 
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and 


ir~ 


13(0,0) 



3 CZ-^-m^)^ 

8 ^4 


1 

128 Z^ 


{zV^ + 2ZW 


+ 3m‘+)6^ 


(A30a) 


In a similar way, it is found that 


13(2,0) « - ^ (Z^-m^)62 (A30b) 

l3(-2,0) « - (Z‘+-m^)62 (A30 c) 

13 ( 4 , 0 ) « — ^ (Z^-2Z2m2+m^)62 (A30d) 

^ ^ o ^ 

l3(-4,0) «=; — ^ (Z'*-2Z2ni2+in‘+)62 (A30e) 

256Z^ 


Further, it is seen from equations (A26) and (A25) that ^3 0 “ 0 for any 
values of Z' except Z' = Z, Z ± 2 , and Z ±4. 


Case (b) : 




m +2 m 
^Z' Z 



P^;2[x2 


(1 


x 2 )]P™ dx 


The integration of equation (A27) over (-1,1) reduces the problem to the 
integration of equation (A25) over (-1,1) with Z' = Z, Z ± 2 , and Z ±4. The 
results are 


13 ( 0 , 2 ) . 


^ , (Z‘* - m^)s 2 
64Z^+ ^ 

(A31a) 

13 ( 2 , 2 ) « — i-T- (Z + m)2(Z2 - Zm + m2)62 
961 "^ 

(A31b) 

1 3 (-2,2) «« - (Z - m)2 (Z + Zm + m2)62 

(A31c) 

13 ( 4 , 2 ) « - (Z+ m)3(Z- m)62 

(A31d) 
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(A31e) 


l3(-4,2) « (Z-+ m) (I - m)362 

384^4 


Case (c): I^^g(AZ,-2) = 


/.! 


^[x2(l - x^)]P^ dx 


To integrate over (-1,1) equation (A28) reduces the evaluation of 
g(AZ.,-2) to repeated integrations of equation (A25) over (-1,1). Again, 

only the results will be given: 


l3(0,-2) « 

13(2, -2) « - (Z-m)2(Z2+^„,+n,2)62 

l3(-2,-2) ^-r- Cl+m)^(t^-lm+m^)6^ 

961“^ 

13(4, -2) « (^-m)3(Z+m)62 

l3(-4,-2) (l+m)3(2-m)62 


(A32a) 

(A32b) 

(A32c) 

(A32d) 

(A32e) 


Case (d): ^(A^,-4) = f P®;‘*(l-x2) 

’ •'-1 


2^2pm ^ 


These integrals can be evaluated by integrating equation (A29) over 
(-1,1). This reduces the problem to repeated integrations, over (-1,1), of 

which is obtained by dividing equation (A28) by x^. The expressions 
obtained are : 


13 ( 0 , -4) 


13 ( 2 , -4) 


256^4 


(Z-m)2(Z+m)262 


384^4 


(Z-m) 3 (Z+m)62 


(A33a) 

(A33b) 
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l3(-2,4) 
13(4, -4) 
l3(-4,-4) 


384^4 


(Z-m) ^ (Z+m)6^ 


6x256^4 


4^2 


6x256^4 


(Z-m) ^6 


4*2 


Cl+m)^S 


Case (e): I3 q(AZ,4) = f - x^) 2p^ dx 


(A33c) 

CA33d) 

(A33e) 


Evaluation o£ these integrals by methods shown earlier gives 


13(0,4) 

13(2,4) 

l3(-2,4) 

13(4,4) 

l3(-4,4) 


1 


(Z^-m^) 2^2 


256Z4 

^ (Z+m)3(Z-m)62 


384Z4 
1 


384Z4 

1 


6x256Z4 


6x25614 


(Z-m)3(Z+m)62 
- (Z+m)462 


(Z-m)462 


(A34a) 

(A34b) 

(A34c) 

(A34d) 

(A34e) 


This completes the evaluation of I3 for all allowed combinations of Am and 
AZ, and makes it possible to write expressions for the overlap integral 
(eq. (Al)) 


6(v,Z,m;v+Av,Z+AZ,m+Am) = ginh ' " S 2Lv^^ (I2+I3) 


For example. 


6 2 ^ v+ 1 ^ 

6(v,Z,m;v+Av,Z,m) = B(Av,0,0) = ' 


C3Z4+2Z2m2+3m4)62 

4Z2 128Z4 
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COMBINATION OF THE (Av, A^, Am) 


The expression for the overlap integral is dependent on m (Z,. . .,-2,-1, 
0 = m = 1,2,. . .,Z) and yet the energy of the molecule is independent of m, 
which only describes the orientation of the angular momentum vector. There- 
fore, an average, over all possible values of m is more appropriate for use 
in energy exchange problems. However, it is not 3 but 3^ that is propor- 
tional to the transition probability, and the average 3^ over all possible 
m will be sought. 

This average 3^ will be found subject to the assumption that Z >> 1, 
since l 2 (AZ,Am) and l 3 (AZ,Am) have been derived on the basis of this 
assiimption. 

3^(Av,AZ,Am) = 2 ~ zVl 3^(Av,AZ,Am) (A35) 

m=-Z 

The sum on the right of this equation can be found exactly, but the same 
results are arrived at by writing 


3^ (Av, AZ, Am) 


1 

2l 



3^ (Av,AZ,Am) dm 


The average 3^ is, of course, independent of m. 

Instead of deriving 3^(Av,AZ,Am) for each possible combination of AZ 
and Am, one such derivation will be made, and the result of this same proce- 
dure applied to the other combinations will be stated. The example chosen is 
3(Av,0,0), whose square is 


b2(av, 0,0) = K ( — 1 (•3J4+2^2In2+3m‘*)]62 

\ 161^ 256J^ 

+ — ^ C9Z®+12^6m2+i8Z‘+m‘*+12Z^m6+9mS)6'+) (A36) 

128^8 ) 


where 


Then 
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M.. 


lllll■lll■ll lll■l■■llll 


I im IIIIIIIII Min ■ ■ 


b2(Av,0,0) 



(Av,0,0)dm 



11 

8960 



(A36) 


In a completely analogous way, the averaged squares below can be found. 


g2(Av,0,2) = g2(Av,0,-2) = B^(Av,2,0) = g2(Av,-2,0) = K 


U^-3'5 2^-3-5-7 


27 . 32 . 5 / 
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g2(Av,0,4) = g2(Av,0,-4) = b2(Av, 4,0) = g2(Av,-4,0) = K 


V 29 . 32 . 5 . 7 / 


(A38) 


g2(Av,2,2) = = K 


/ 31 ^ 36762 ^ 56336 ^ \ 

\28. 5 ^ 29.3.5.7'^ 2 12. 34. 5.7 j 


(A39) 


g2(Av,2,-2) = b2|;av,-2,2) 


(- 1 —+ 1762 ^ 3836^ \ 

\ 2 ®- 5 ^ 2 ^- 3 - 5 . 7 ^ 212.34.5.7^ 


(A40) 


g2(Av,2,4) = g2(;Av,4,2) = e^(Av,-2,-4) 


62(Av,-4,-2) = K 


/ 2336^ \ 

\2l5.3‘^.7/ 


(A41) 


g2(Av, 2,-4)= g2(^v,-4,2) =F(Av,4,-2) =^(Av,-2,4) = K 
F(Av,4,-4)= ^(Av,-4,4) = K 
?2(av,4,4) =F(Av,-4,-4) = K 


(A42) 

(A43) 

CA44) 


Each of the squares g2(Av,A7,Am) is proportional to the transition 
probability P(Av,AZ,Am). The total transition probability, for all possible 
values of Am is given by 


P(Av,A 7) = P (Ay, AZ, Am) 

Am 

Consequently, the constant, g2(^Av,A7) to be used in computing this transition 
probability is 
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where 

_ ( \^ rnftoj(v+l,0)l ^ ( 6^ \^/ v+l,0 \ 

Vsinh 6/ L 16uE(, J \sinh 6/ \4L2a / 

and (v+1,0) = v+1 or v accordingly as Av = ±1. 
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APPENDIX B 


CALCULATION OF COLLISION IMPULSE FOURIER TRANSFORMS 


A slightly asyiranetric collision perturbation of the form 

U(t) = Uoe^^ sech^ at (Bl) 

has been used to account for the effects of irreversible excitation of 
vibrational energy during collision. The dimensionless Fourier transform of 
this impulse is 


/ oo 

e^^ sech^ (at) e^^^ dt 

00 

This integral, which can be written 

/ oo 

sech2 X dx 

-OO 

will be evaluated by contour integration. Accordingly, we integrate the 
complex integral 


(B2) 


(B3) 


sech2 z dz 

about the rectangular path shown in the sketch. The integrand is analytic at 
all points interior to this contour except at the pole z = (7ri/2) . 
Consequently, 


(-C+ 7Ti) 


iy 


(C + iri) 


irt 

2 


-C 


0 


c 


X 
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f sech2 x dxtif* sech^ (ciy)e';"“/»> dy 

•'-c •'o 

+ r"" sech2(x + 

•'c 

+ ij* ^ ^^^^^sech^(-c + ^ c+iy) jjlesidue at (B4) 

When we take lira of both sides of this equation, we get 

C-^oo 

/ oo 

sech2 X dx [i _ e 


+ i lim 

C^QO 


+ i lim 

C->C)0 


r sech^Cc + iy)eCWa)(c+iy) 

•^n 


dy 


r 


(q/a) C-c+iy) 


sech2(-c+ iy)e^^“^^^^ ^ c+iy) jjy _ 2Tri jResidue at 


(B5) 


Since sech^ z = 1 + tanh^ z, and lim tanh z = 1, it is clear that 

z~yoo 


lim 

C“^ oo 


► 0 


J gCq/a)C-cH-iy) 3ech2(-c-My)ef"<"/^^f-'^"i>^^ 


dy 


^ lim TT 

C-^00 


^(q/a) C-c+iy) + tanh^C-c + iy) jg 


0 (B6) 


The second integral in equation (B5) is treated by writing its integrand in 
terms of exponentials. 




4 g [ (q/a) - 2 + (iw/a) ] (c+iy) 
1 + 


(B7) 


59 



Since (q/a) << 1, it follows that Cq/a) - 2 < 0 and therefore 
lim g(y) = 0. Hence, 

C-H» 


lim 

C->oo 


/' 

•'n 


^(q/a)Cc«y) * iy)^CWa)(c*iy) 


<. lim ir|g(y) | = 0 (B8) 


and equation (B5) becomes 


>% oo 

f sech2 X eti”''®)* dx = 

-.oo 


i ^Resi 


2Tri|Residue at ~ 


[^1 + e(7T/a)(qi-o)) 


(B9) 


In order to find the residue of 


fCz) = sedi2 2 


at z = Tri, it is of help to write 


fCz) 


e[(q-i(o)/a]z ^ 

cosh^ z q(z) 


where both p(z) and q(z) are analytic at z = iTi/2, and can therefore be 
expanded about that point. 


f(z) 




'q"C^i/2) q'" (7ri/2) 

fz - 


2 6 

\ 2) ■ 



If we carry out the division of the series, we obtain 


(BIO) 


P (~) ■ 

2p* (iri/2) 

2 p(Tri/2)q'” (7ri/2) 1 

^ \ 2 J (. 

r'(ni/2) 

3 [q"CTTi/2)]2 1 



/_ TTi \ 

V" ■ —) 


- -) 


(Bll) 


and the residue is given by 
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Residue at = 


3 


(x) 


= + io) ^ ^ r Cq+it»/a) 1 CTri/2] 


Finally, 


^ -2a)iri / q ■*• ia) \ 

p. , « f sech2 X dx - ° N -?> 


-2g)7ri ( °i * \ ^[(q+iw)/a] (TTi/2) 


[ (q+iaj)/a]iri 


(B12) 


arid for (ucD/a) >> 1, this gives 


p* 


27TW 


((0 - iq)e-^^“/2a) 


SO that 


|p *|2 ^ 4^2 g (Tra)/a) ^ 

The factor [1 + (q^/oo^)] represents the correction to the usual adiabatic 
collision transition probability, which is due to the asymmetrical factor e*^^ 
in the impulse function. 
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APPENDIX C 


VIBRATIONAL ENERGY AND ITS DERIVATIVES AT THE 
COLLISION TURNING POINT 


Consider a diatomic molecule subject to the perturbation 


U = Ug + j (grad U)g y 


(Cl) 


where y = (p - oscillator coordinate. When the molecule is ini- 

tially at rest (y^^ = 0) > the maximum amplitude of classical motion induced at 
any time t is, according to Green’s method 


m 




(grad U)e 

2yo) 


cos ojT dr + 




(grad U)e . , 

PS Sin (jJT dT 

2yo) 


so that 


m 


/ 


^ U)e icot 

2yu ® 


dt 


(C2) 


The net amount of vibrational energy excited is, in units of hto 

I -1- 

1 


e = 


AE 


hw 2hco 


Syhoi 


r (grad U)g dt 

•'-CO 


(C3) 


This classical result is exactly the same as the quantum probability of 
transition 0-^1 given by small perturbation theory 


e = P 


01 = 


1 f* „ llgt 


dt 


(C4) 


where Hqi is the harmonic oscillator matrix element 


Ho. = (grad 


U). 


(C5) 
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When the transition probability becomes too large for small perturbation 
methods to apply, Kemer's solutions of Schrodinger * s equation are used. The 
probability of transition m ^ n due to a forcing function acting on a 
harmonic oscillator is found to be (ref. 24) 


P^(t) = m!n! 


(C6) 


where the polynomial is 


^mn(^) “ 


m<n 

z 

3 = 0 


(-l)^e’^. 


(n - j) ! j ! (m - j) ! 


(C7) 


For present purposes, transitions from the ground state are of interest. 


n -e 

> = ^ ^ ■ 
on n! 


(C8) 


Note that Pqi reduces to the usual small perturbation result (eq. (C4)) 
when e is small. Also the sum of all is unity as it should be. The 

total energy in all vibrational modes is 


03 



n=o 


CO 


E 

n=l 



hwe 


(C9) 


Where the interaction potential is exponential, (grad U) = U/L and 
equation (G9) can be expressed 

B rt . ^ 

^ Ue(t)e^“^ dt 


V(t) ^ ^ _ 
ho) 


CCIO) 


which is the result used in equation (109). Hereafter, we shall drop the 
subscripts on 3 and U for convenience. 


The vibrational energy at t = 0 is 


/ 


2 


2 ' 

Vo = -^1 
° hi 

I U cos (ot dt 

+ 

1 U sin (lit dt 


1 

CO 


— oo 



(Cll) 
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The first integral on the right is negligibly small compared with the second. 
The latter can be treated by repeated integration by parts to give 



U sin tot 


dt 




Ur 



+ . 


0) 


(C12) 


where we have used the fact that U and all its derivatives vanish at 
|t| = ". Thus, to terms of first order in to“^ 


Vo = 


hto 


(C13) 


The derivatives of V at t = 0 are obtained by differentiation of 
equation (CIO). The cosine transform of U and all odd derivatives of U 
at t = 0 are taken to vanish 


/ 


0 

oo 


U cos (ot dt = Uq 



(C14) 


which is equivalent to neglect of asymmetrical terms in U and V. Then the 
results are 


Vo = Vo = Vo = 0 


"v - U ii 

''O ^0^0 


(CIS) 

(C16) 


The vibrational energy at the turning point is seen to vary essentially 
as t^*. These results are valid for any impulse function U(t) where the 
duration of impulse is long compared with provided only that U varies 

exponentially with distance between the collision partners. Similar relations 
obtain in the more general case, only U is then replaced by (L grad U) . 
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